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A THEORY FOR DETERMINING OPTIMUM DYNAMIC SYSTEMS 


N. I, Andreev 


(Khar’ kov) 


Sufficient conditions for an extreme of a certain functional, taken as a criterion 
for comparing various dynamic systems are derived, 


Two examples illustrating the proposed technique for-choosing an optimum 
dynamic system are given, 


It was shown in [1] that the problem of finding an optimum dynamic linear system that would ensure an 
extreme in functional I, which represents a given function n+ 1 of quadratic* functionals of a unit pulse trans- 
fer function of system k(t), can be reduced to finding an extreme of the quadratic functional 


Tyy = 6,1, ()] +... + n + lk 


Function kg (t) which ensures an extreme in functional Ip, depends on the artificially introduced co- 
efficients @;. The values of coefficients 6; = 6; 9 at which the extreme of functional I[k(t)] is obtained can 
be determined by one of the two methods suggested in [1], 


The required conditions which must be satisfied by k’(t) in order to obtain an extreme in functional I 
consist of the following: 


1) = Ko (t, O49, » Ong) Must be the extreme, i,e,, must be the solution of Euler's equation corre- 
sponding to functional Ip; 


2) al 


or (a special case) 


a, 7% 


In [1] the question of sufficient conditions which function k’ (t) must satisfy in order to ensure an extreme 
in functional I was not raised, 


* It should be remembered that the term quadratic functional is applied to a functional ly when it possesses 
the following properties: 
1, (0) + AK = Ay + AA + 


where k(t) and x (t) are arbitrary functions from a class of permissible functions, and Ag, Ay, and Ag are 
numbers related to functions k(t) and «x (t) but independent of A, 
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The present article gives a generalization of the material presented in [1] and in addition the sufficient 
conditions for an extreme in functional I = I[k(t)]. 


Let us formulate the posing of the problem, 
Let us assume that functional I[k(t)] represents a given differentiable function © of functionals I ,{k (t)], 
++ Ingg(k(t)) (these functionals need not be quadratic): 
IT=@(I,,..., In4s). (i) 


Let us also assume that there is a method of determining a function which would ensure an extreme in 
functional 
with any values of coefficient 6; and that in particular the required and sufficient conditions which must be 
satisfied by function ky (t, 6y,..., @n) for ensuring an extreme in functional ly are known, 


With these conditions it becomes necessary to state a method of determining function k’ (t) which ensures 
an extreme in functional I, and to find in particular the required and sufficient conditions which function k’ (t) 
must satisfy, 


In order that function k’ (t) should ensure an extreme in functional I it is necessary that the first varia- 
tion of this functional with k(t) = k(t) should equal zero: 


dl (t) + Ax 
= = Q, (3) 
where « (t) is any arbitrary function from a class of permissible functions, 
Taking into account (1) let us rewrite condition (3) in the following form: 
OD dl, (k(t) + Ax (4) 


The first variation of 1, as it will be seen from expression (4), can equal zero if at least one of the 
following conditions is fulfilled: 


1) 6M /0I;=0, iqi,...,n+1; (5) 


2) if at least one derivative 3 / 01; does not equal zero (without reducing the generality let us 
assume that 04/@1p4, # 0) and if the following equation holds 


dl, [k° (t) + AxK(t)] dl, (t) + AK(t)] 
“ dd 4-0 


+ (t) + AK ‘t)] (6) 


/ 
= 
/ 


Equation (5) is, of course, a necessary condition for an extreme in function $(I,,...,In4,) with 
certain values of arguments 1; =1j 9. 


If with a certain function k (t) = kg (t) the functional I; assumes the values lj =1j9, it should be checked 
further whether the many variables function @ satisfies the sufficient condition for an extreme with values 
I; =1j9. This sufficient condition is the inequality [2] 


where 
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51,51; <0 


with arbitrary values of 61; and 61; variations, 


If conditions (5) and (7) are fulfilled for a certain function kgg(t), this function will ensure the extreme 
for functional I, 


In problems arising from practical requirements condition (5) is not usually satisfied for the selected 
(given) class of function k(t), Functional I becomes in such problems the criterion for comparing various 
systems (for instance systems of automatic control) and represents a smoothly changing function of functional 
I,, moreover the extremes of function (1) are usually outside region G of the variations of functionals l; = 


=1(k(t)). 


We shall consider henceforth that in the region of variations of quantity 1; [for a given class of function 
k(t)] all partial derivatives 66/ 01; cannot equal zero simultaneously, In this case, as it has already been 
pointed out, the required condition for an extreme in functional I is condition (6), 


Equation (6) is a complex integro- differential equations we shall solve it, therefore, in two stages, First 
we shall find a general form for the solution 


df, [ky (t) + dT, [ko (t) + 
[ko(t) +4 «(0)] (6a) 
+ ™ 0, 


where are as yet undetermined parameters, This solution kg(t) depends on parameters : kp = kg (t, 
6»). 
Next let us determine 6; = @j9 in a way to satisfy condition 


/ aI, (8) 


= 4; i, d= i,... ne 
OD / 91,45 (helt, » Op)) 


Such a splitting into two parts of the solution of equation (6) is possible because the solution of this 
complex equation coincides with the solution of equation (6a) if the coefficients of (6) are equal to 


a,=%, (9) 


Hence, the solution of k’ (t) lies in the class of functions kg(t, @3,..., On» The values of coefficients 
6; with which condition (9) is satisfied are found from condition (8), 


Equation (6a) is the necessary condition to obtain with k(t) = ky (t) the extreme of functional ly [see (2)], 
According to the second condition of posing the problem the method of determining the general solution of 
equation (6a) is known, If this solution k(t, 6), ..., 0) of equation (6a) can be determined in its final 
form, we shall obtain a system of algebraic or transcendental equations, determining the value of @j9, by 
substituting expression kg(t, 63, ... ,» @,) im the left-hand sides of equations (8) or, which is the same thing, 
in equation (9), Usually under these conditions the solution of the system of equations (8) is not single- valued, 
In such cases the value of functional I should be calculated for coefficients 6; which correspond to the varia- 
tions in the solution of the system of equations (8) and those in the neighborhood of the variations, selecting 
from them the variation at which functional I has the largest (smallest) value, 
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In practice there can be cases when the solution kg(t, 6;,..., 9p) in a certain region of variations of 
coefficients @; has one analytical expression and in another region of variations a different one, In such a case 
the solution of the system of equations (8) should be sought for each of these regions, Moreover, in one of the 
regions there may not exist an extreme for functional Ij, Such a region of @; parameters should be excluded 
from further consideration, 


In searching for coefficients @;, which satisfy the system of equations (8) even greater difficulties may 
arise, Sometimes, it is more convenient to determine coefficients by another method, This method consists 


in substituting solution kg (t, 63, . . 9p) of equation (6a) into functionall, In this case functional I is con- 
verted into a function of n variables 6;,...,9n2 I=F(0,,..., 4). 


Values of coefficier.ts 6{9 which correspond to the extreme of the functional I are determined in this 
case by the system of equations 


OF (01,....0,) 
00; 


(10) 


Moreover, one or several of the equations (10) can be substituted by the condition (conditions) that 


+» On) does not exist, 
004 


The cases when OF / 06; do not exist are considered to be exceptional and will not be dealt with any 
further, 


If a combination of parameters of 6; 9 corresponds to the extreme of 1 =I [ky (t, 04,... , Op)bit is the 


solution of the system of equations (8) and (10), since these systems of equations are the necessary conditions 
for an extreme in functional I, 


If the combination of parameters 6; 9» however, does not correspond to an extreme I =1 [kg (t, 6;,..., 
=F Op) but is the solution of the system of equations (8), it coincides with the solution of the 
system of equations (10) and,vice versa, the solution of the system of equations (10) coincides with the solution 
of the system of equations (8) with the condition that 


In order to prove this statement let us assume that condition (8) for function kg (t, 64, ..., Op) has 
been fulfilled with a certain combination of the values of parameters 6; =09j9, i=1,...,m. Let us show 
that with the same combination of parameter values 6; =@;9 condition (10) is being fulfilled for the same 
function kg (t, 63, . . «.. 9p) which satisfies the necessary condition for an extreme in functional lj, This 
necessary condition as it was pointed out before, is the equality to zero of the first variation of the functional 


= % ol (t)] + + Onoln (t)] + Inga fh (0)] 
when k (t) = kg (t, Oy. Ong)- 


Hence, the conclusion that the partial derivatives of ly with respect to parameters 6; with 6; = 
= jo (i, j=1,..., 0) are equal to zero 


a, ALA 
i 950 i 95 =0j9 
"70 


00; a0; 
- 
00, 00, 
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By substituting in relationships (12) coefficients 649 in front of the partial derivatives 91, / 06, by the ratios 
of the partial derivatives (8) and then multiplying relationships (12) by 8¢/8In4, we obtain equations (10): 


aD am OF 
13 


where k(t) = kp (t, 


Let us show that it holds conversely as well, if condition (10) is fulfilled by a certain function kg (t, 64, 
. «+» 9p) which is the solution of equation (6a) with parameter values 6; = @j9, then these parameter values 
will turn equations (8) into identities, 


Equation (10) can be written in an expanded form (13) or divided by the derivative 04/981,,, and 
written as 


ad / al, Ay, Ay 


om 0, (14) 


0; = Oj, i, 


By considering that in this case equations (12) still hold, and comparing the systems of equations (12) and 
(14) it is possible to conclude [3] that if condition (11) is fulfilled the values of 6; =@49 also satisfy equations 
(8), In practical problems condition (11) is normally fulfilled, Hence, conditions (8) and (10) determine identi- 
cal combinations of values 6j9,i=1,...,0. 


For future reasoning it should be noted that equation (6a), which is the necessary condition for function 
k’ (t), corresponding to an extreme for functional I, can be obtained from two considerations, Functional I of 
function k(t) is also a function of n+ 1 variables Ij, i,e., is a function of a point in space of n+ 1 dimensions, 


Region G of the n+ 1 dimensional space, in which function  (,. .. , In4q) is given, is determined 
by the limits in which the values of functionals 1j[k(t)] vary with a given class of permissible functions k(t), 
If function #(l,...,In+y) reaches an extreme within the region of variations of variables lj, the values of 
variables 1; 9 at which the extreme is reached are determined by conditions (5) and (7), 


Other extremes of functional I corresponding to the larges, or the smallest values of function ¢ ee 
In+1) can be obtained only on the boundaries of region G, Hence, in order to find these extremes of functional 
I it is necessary to find the boundary of region G and then to find the largest (smallest) value of function @ on 
the boundary, 


The boundary of region G can be determined the following way, Let us assume certain values for functionals 
In: 


With these conditions let us determine the extremes of functional 1,4; =1p4,[k(]. It will be seen from 
[4] that in this case the problem of determining function k(t) which would ensure a conditional extreme for 
functional I, is reduced to finding a function which would ensure an absolute extreme for functional 


I, = 6), + + + 


where 6; are the undetermined parameters, whose values are found from conditions (15), The necessary con- 
dition for function kg (t, 64, ... ». @p) to ensure an extreme for functional ly is equation (6a), 


It should be noted that in this posing of the problem the values of 6; , which satisfy condition (15), are of 
no interest to us, It is important for us to know that functional I; attains its boundary values when k(t) = ky (t, 
61,...,0,) where function ky is the solution of equation (6a), Whence, we conclude that the boundary of 
region G may be determined by the relationship 


= [ko (t, p41 (16) 


| 
a 
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and represents an n-dimensional space, Functional I represents on this boundary a function of n variables 6;: 


The necessary condition for an extreme of function F is written in the form of relationships (10), 


The above new approach to the determination of function ky (t, 63, . .. , @p) and the values of coefficients 
6; which correspond to the extreme of functional I are not related to the differentiability of function # (, 
«+ + » 1,41) with respect to values of I;, This approach to the problem can also be used when selecting an 
optimum nonlinear dynamic system which ensures an extreme for function @ = (ky,..., ky,) where k; are 
values related to the characteristics of a nonlinear dynamic system, The suggested approach to the solution of 
the problem of finding an extreme for functional I is more convenient when sufficient conditions for the ex- 
treme of this functional are found, In this connection let us examine two cases, 


In the first case, when theextreme of (1,,..., 143) is reached inside region G of the variation of 
functionals I;, the sufficient conditions for the extreme are represented by relationships (5) and (7). 


In the second case let us assume in addition that @ (I;,... , 1,44) is a monotonically decreasing 
function of I, in region G and the boundary of region G is a convex surface, and then split the problem into 
two s’%ges, In the first stage we shall obtain sifficient conditions for ensuring that with k(t) =kg(t, 63,..., 
where [6; — < 4 being positive numbers, points (;,..., 1,43) lie on the region G boundary, 
In the second stage let us obtain sufficient conditions for function F (@;,..., On) =I[kp(6y,...,6p)) to 
reach the extreme on the boundary of region G with certain values of parameters 6; = 6; 9. 


With k(t) = kg (t, 63, ..., On) where (6; — Ojo] < points , lie on the region G boundary 
if functional I,,, reaches the extreme under conditions that =Cy wherei=1,...,mjand dy < di, 
di, and dj, being numbers determining the region of number Cj variations, The numbers dj, and dj, are 
functions of and Hence, in order to ascertain that with k(t) = kg (t, 6;,..., 0p) points (y,..., 
I,41)) Me on region G boundary it is necessary to have sufficient conditions for functional Iy to reach the ex- 
treme with k(t) = kg (t, 0,....6,) where [6; — 6j9] < €j. When formulating the posing of the problem 
at the beginning of the paper, it was assumed that these sufficient conditions for the extreme of functional I; 
are known, 


Now the problem of the second stage must be solved, i.e., the sufficient conditions for function F (6;, 
+ e+» On) = Ike) to reach the extreme at 6; =6j9 must be found. Such conditions have already been written 
for function © (y,...,1,41). These conditions (5) and (7) assume the following form for function F 


OF 


(17a) 
6,08, 66, 50; > 0 


i, j=1 
(17b) 


0, 00, 90; <0, 


i, j=1 (17c) 


where the partial derivatives are calculated with parameter values 6¢ =@e9, €=1,..., m, Which satisfy 
conditions (17a), and one of the conditions (17b) must be fulfilled with any values of parameter 6 6; and 
66; variations, 


Thus, the sufficient conditions for an extreme of the functional consist of sufficient conditions for an 
extreme of functional 1, andconditions (17a) and (17b), The more complex conditions are usually those for 
the extreme of functional 1,, These conditions depend on the form of functional l;. 


Let us now give an example often encountered in practice when selecting optimum linear systems of 
automatic control, Let functional I; be quadratic of the form 


or 
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(18) 


where R(t) is the correlation function of a stationary continuous random process, k(t) the integrated function 
in the interval 0 — T, and Cj, Dj are some numbers, 


The necessary condition for an extreme (in this case a minimum) for functional ly can be obtained by 
means of making the first variation of lL, equal to zero, In this concrete example it is easy to see that this 
condition has the form: 


m T 
—t) k(t) dz i \ de = 0. (19) 
Re )k(x)d +> [e ik (x) d + 0 19 


Let us show that this necessary condition for the minimum of functional Ip is also a sufficient condition 
ifC; = 0 wherei=1,...,m, For this purpose let us calculate the difference 


Ip (t) + (t)] — [9 (2)], 


where k°(t) is a function satisfying condition (19), i.e,, an extreme, « (t) an arbitrary function integrated over 
the time interval 0 — T, 


It is easy to see that 


T 


+3 th de + dh + 


0 


i—o 


since the first term is a variance of a random process, and the second a sum of quantities which are not negative, 
Thus, the sufficiency of conditions (19) have been proved, If some-of the coefficients Cj < 0, however, con- 
dition (19) becomes insufficient in a general case, In practice the coefficients of Cj are usually positive, 


The above method is also applicable when function k has two variables k = k(t, 1) and I is an operator: 
1=I[k(t, r)] =1(t), In this case variable t should be considered a parameter, In this paper as in [1] for 
simplifying calculations parameter t has been omitted, 


In conclusion, we give two examples which illustrate the technique and basic principles of this paper 
and of [1]. 


Example 1, Let a signal y, consisting of a useful signal » and interference z, be impressed on the input 
of a linear automatic control system. The useful signal x can be represented in the form 


X + vt, 
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where xX, is a random ‘quantity with a large variance, v is a given quantity and t is time, 


The interference is a stationary random process subject to the normal distribution law, with the following 
probability characteristics 


mathematical expectation =0, 
correlation function R,(r) = 
Quantities x and z are not interrelated, i,e,, M[xz] =0, 
, The process at the output of the automatic control system is represented by expression 


=| 


where 6 (t) is the law of change of the system output coordinate; k(t) the unit pulse transfer function of this 
system with 0 <t <T, k(t)=0 when t< 0 andt> T, 


Let us assume that function k(t) is subject to the following limitations, 
k(t) is the integrated function 


T 
k(x) dt =1. (21) 


It is required to determine the unit pulse transfer function k(t) at which the maximum probability is 
reached that the error of extrapolation 


= 2+ %)—| 


will not exceed in its modulus a certain small quantity «, Here Te, is a given extrapolation time, 


The probability that | A| < €, where ¢ is a small quantity, can be determined with great accuracy by 
the following approximate formula: 


e (23) 


2e 
V 


where a=M [A] is the mathematical expectation of the error of extrapolation, o? = M[(A— a)*] is the 
variance of the extrapolation error, 


Henceforth, we shall not operate with the probability P,, but with quantity P, proportional to P,: 


P= 


Vin 
Qe 


(23a) 


It will be easily seen that by changing from P, to P no substantial change has been introduced, 
Quantities a? and o” are functionals of function k (t) 


a? — {M[A}}? = E (t + Te) y(t — t) k(x) az [re + ( tk (1) (24) 
0 0 


= — =m | = ( R, (A — k (a) k(t) dd de. 
0 0 


0 
9 
| 
| 
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In this example functional I corresponds to quantity P, functional I, to functional a®, functional I, to 
functional o”, Functional I; is in this case a quadratic functional, since a* and o? are quadratic functionals, 
Hence, in this case functional I; has the form: 


= 00? (26) 
Taking into account (24) - (26) we can write: 
TT T 
a) + [7, th (2) ax’. 
00 0 


Function k(t) = ky (t, @) which ensures an extreme (minimum) for functional l;, with @ > 0 and con- 
dition (21), ensures an absolute extreme for functional 


00 0 0 


where y» is an undetermined factor found from condition (21), 


A necessary and sufficient condition for function kg (t, @) which ensures a minimum for functional I is 
the integral equation [5, 6] 


T r 
6 \ R, (t — t)k (+) dt + v% \ tk (t) dt +v°T, t—y, = 0, (28) 
0 0 


which is a particular case of e mation (6a), 
The solution of this equation from [7] has with R,(r)= e" It! the form: 


ky (t, 8) = Ay + Ayt + C,8 (t) + (t —T), (29) 
where § stands for a delta function, 
Coefficients A», Ay, Cy, and D, are determined by a system of equations obtained after substituting (29) 
in (28): 


T 
+ Are + C18 (2) + + 
0 
T 
t[Ay + Ait + (x) + DB —T)] dt + — = 0. 
0 


After integration, having equated to zero coefficients of fot, eM, eM and using condition (21) we 
obtain five equations for the determination of five unknownsA , Ay, Cy, Dy, and yo: 


a Apt(eT+1) A — aD, =0, 
a A,— A, —a@C, = 0, 
T At A+ Di=0, 
T T? 20 
+) A, + TD, =—T,, (30) 
20 
Ao — ¥o= 9, 
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(ar +1) 7 + a|+(ar +S) #7, | (31) 
9 
a? (ar + +09 (aT + 2)T, 
A\=-- De 


where D, is a determinant consisting of the coefficients of the first four equations of system (30) and taken 
with a reversed sign, i.e., 


Dy = 40S +207 +1) 4 20878 (32) 


T 
We shall now require to express the integral vy = \ tk(t)d+ in terms of the initial values of the problem, 
0 


Omitting intermediate calculations we obtain 
T 


v= tk, (t, = — A, —T,, (33) 
0 


In order to fulfill the second stage of the problem it is necessary to express a” and o” in terms of the 
initial values and coefficient 
By using formulas (24), (31) - (33), we obtain 


a? (0) = + v= 41) = 


=4%_@ 


[os + (055 +1) +2007 4 2 are 4 


Using formulas (25), (28), (31) - out we obtain 


Dy 
2a? [7 (14-27) 4 (ar + —aT (1447) a(ar +2)7 


= 2a 


Introducing notations x= «a7T,kK,=«7,, y?= Ss , after transformation, we obtain 


+ K(2 + 2x+ + +0 +K)Ke(Ket (35) 
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Hence 
A 
: 


With these notation expression (34) assumes the form 


(n+ 
a? (0) Dy 6?, 


where 


2, 2 x4 
Do = 2° (2 + K) + 2k + 2K + 
In subsequent calculations, we assume 
kK=3, = 20. (36) 


With these concrete numerical values of the initial values the quantities a? and o* can be expressed 
in terms of parameter @ as follows: 


Bt 2.8 6? 0.40 6? + 0.385 0 + 0.65 


According to the proposed method it is now necessary to find the value of parameter 6 = 6», which 
solves equation (8), In this case (8) has the form: 


OP | da* 
or 
20 


Substituting in equation (38) expressions (37), we obtain the relationship 


2,8 0? 
0,40 6? + 0,385 6 + 0,65 


i—6, 
which after transformation has the form: 


6 + 76? + 0.66 6 — 1.63 = 0. (39) 


From Descartes theorem [3] it follows that equation (39) has one positive root, The other two roots are 
negative, These roots are not taken into account since it can be shown that with negative values of parameter 
@ the sufficient condition for a minimum in functional ky [see (20)] is not fulfilled, and that a maximum in 
functional I = P is not reached either, 


Before attempting to find the positive values of root 69, of equation (39), let us note that this root is 
smaller than one,since-it is only with @ < 1 that the mathematical expectation a(6) decreases and a rise in 
functional P can be obtained. Hence, the root of equation (39) should be sought in the interval 0 to 1, The 
value of this root is equal to 6, = 0,42 with an accuracy of 0,01. 


We do not propose to check the conditions for a maximum in functional P, instead let us determine the 
maximum for P = P(@), by plotting the graph of this function, 


By using formulas (23) and (37) a®, o, P, and P* were calculated for several values of @ and the table 
given below was compiled, 


In addition to the values of P(@) and P%(@) the table also shows the ratios between these values and P (1) 
and P*(1) which correspond to the optimum system of automatic control, determined by the criterion of the 
minimum mean-square value of the error, 

a* 


Functional P* has a completely definite physical meaning, By means of functional < Pt = Ss 
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#0) (38) 
| 
a 
| 


it is possible to calculate with great accuracy the probability that |A,| <« and | A,| < €, where ¢€ isa 
small quantity , and A; and Ag are errors in two similar automatic control systems which displace the final 
control element in mutually perpendicular directions, when independent from each other signals with similar 
probability characteristics are impressed on the system inputs, 


6 1.00 0.80 0.60 0.50 0.42 0.35 0.20 0.00 
P*(0) | 0.155 | 0.210 | 0.258 | 0.285 | 0.295 | 0.290 | 0,230 | 0.100 
P*(0)/ P2(1)| 1.00 | 4.35 | 1.66 | 1.83 | 1.90 | 4.86 1.48 | 0.65 
P (0) 0.395 | 0.458 | 0.508 | 0.534 | 0.542 | 0.538 | 0.480 | 0.315 
P(6)/P(4)| 1.00 | 1.46 | 1.29 | 1.35 | 1.38 1.36 1.22 | 0.805 


Data graphs showing P(@)/P%1) and P(@)/P(1) were plotted from the table (see drawing), It will be 
seen from the graphs that with ky (t, 6), where @ =@9 = 0.42, the functionals P? and P do attain their maxima, 


A study of the table and the graphs leads to the conclusion 


that the selection of an optimum system in this example by 

a pte means of the maximum P (or P”) criterion leads to higher results 
° N\ wai than the selection by means of the minimum mean-square error 
6 7 SJ criterion, The former gives a P 38% greater than the latter, 
The increase in P* is equal to 90%, 
t2r— 7 “A | It is possible to evaluate functional P* for an optimum 

a system selected by means of the maximum P” criterion, if the 
a. values of aj and o% for an optimum system selected by the 


minimum mean-square error criterion are known, The values 


a6 
Q2 QI Q7 48 a9 106. 


| 
a 
ote! ote 


1 
% 


Certain relations ai, and > of when al, = a? +03 correspond to the optimum 
system selected by the P* criterion, Let us determine maximum P* under condition that a? and o? can assume 


any value and that a*+0? = a} + of, It is easy to see that with this condition the P* maximum equals 
1 


oe + 


of a? and o correspond to the value P? = P§ 


Hence, the value of the functional P? = P?,, which corresponds to the optimum system selected by the Pp? 


1 1 
criterion is smaller than ab 1.€., = ab 


2 
Hence, if ratio -- is small the difference P?, — P? is also small, 
0 


It is also possible to evaluate P%, from below if a3 > o%, For this purpose let us reason as follows, In 
the transition from P? to Pio» a” decreases and o” increases, The values of a* and o” change continuously, 


In the process of their variation a” and o? attain values of a4, and of, at which quantity P? is at a maxi- 
mum, 


Let us assume that a* remains constant and equal to af and that o” increases, With this condition P* 


reaches a maximum when o”= a2, This P* maximum is equal to » hence = . This means 


that if the difference 


is substantial, the changing of the criterion from the minimum mean- 
ea 


0 0 
square error to the maximum P” gives a substantial gain in P*, 
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For estimating the gain in P? it is sometimes more convenient to use the ratio P3,/P}, 


a2 
0 
% (40) 
e ° at /o8>1. (40a) 


It follows from inequality (40) that with small values of a2 /o3 the ratio Pj,/P} is close to unity, i.e., 
that the changeover to criterion P? does not give a noticeable rise in P*, Inequality (40a) shows that with large 
values of ratio ag/ 0} the changeover to the P* criterion results in a substantial increase of functional P* (or P), 


The above evaluations of ratio P%/P? and the formulas (34) and (35) lead to the following conclusions 
with respect to this example, With fixed finite values of ke and y and with x —> o ratio aj /oj tends 
to zero, hence, under these conditions a changeover to the Fs criterion is not expedient, since it does not pro- 
duce a marked gain in P*, With the fixing of the value of y * and Ye Wand when x~ 0 the ratio 
P3, / P}_ tends to oo, and hence with large x, and small x a changeover to the P* criterion is expedient, 
since it ensures a substantial gain in P* (or P), 


Example 2, The conditions of the problem are the same but with the difference that now the value of v 
is not a given quantity, but a random one distributed according to the law of equal probability in the range 
—v, to +vy. The random quantity v and random function z are not interrelated, i.e., Miyz] =0, In this 
example the probability that | A | < ¢€, where « is small, is calculated from formula 


3 
2e 
P. = Pp A = 2 

(| | < e) V x V 30, e dz, 

where oy is the mean-square deviation of the error component at the system output due to the input error z, 

oO; is the mean-square deviation of the error component at the system output due to the random quantity v. 
Variance of is calculated from formula (35) and the variance 0 from formula (34), Moreover, quantity y 

in formulas (34) and (35) is calculated from formula 

2 2 
= =V3 


v 1 


Similarly to the first example, we shall use instead of P, quantity P which is proportional to it: 


v3> 

-= 41 
0 
and in addition to P we shall also consider quantity P*, 


By taking the same initial data (36) as in the first example and utilizing formulas (34), (35), and (41) it 
is easy to plot a graph of function P%6) and P(@) in the range of @ between 0 and 1, These graphs permit 
one to determine the maximum values of P* and P and corresponding ratios P*(@»)/P*(1) and P(6,)/P(1), 

Below are the values of these ratios: 


P(i) 
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P? (0.43) P (9.43) _ 4 99 


Thus, with the initial conditions discussed above the changeover to the P criterion leads to an increase 
in P* of 18% and in P of 9%, With other initial conditions the gain in P (or P*) will be greater or smaller, 
The conclusions drawn with respect to the increment in P or P* with relation to parameters x and x, deduced 
at the end of the first example hold for this example as well, 


SUMMARY 


The paper gives the necessary and sufficient conditions for an extreme in functional I= 6 fate 
«+++1i41 [k(]} for the case when the necessary and sufficient conditions for the extreme in the functional 


+. Onin are known, 

The results obtained in the paper can be used for selecting optimum dynamic systems, 

Particular examples at the end of the paper show that the use of a more complex criterion (P or P”) leads 
in certain cases to an increase in the quality of the sclected system (i.e,, P or P*), 
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MATHEMATICAL SIMULATION OF DRY FRICTION®* 


G. I. Monastyrshin 


(Moscow) 


Simulation of dry friction is considered in three cases: when starting friction is 
equal to kinetic friction, when it exceeds the latte and when it depends on the duration 
of the contact at rest, Tried out circuits for simulating the movements of solid bodies 
with one degree of freedom and simulating minor oscillations of a gimbal gyroscope 
axis are described, A method of simulating, by means of differentiating amplifiers, 
systems whose differential equations have not been solved with respect to the higher 
order derivatives is given, 


Dynamic systems often contain mechanical devices, In automatic control systems, forinstance, mechani- 
cal devices are usually employed in sensing elements and final control stages, In view of this it often becomes 
necessary to simulate nonlinear processes due to dry friction. 


It is usual to base mathematical simulation on the analysis of dry friction given by the French scientist 
Amonton (1663-1705) [1]. It is assumed that at rest the friction force balances out the effort tangential to the 
sliding surface, Movement will result when the above external effort exceeds the maximum possible value 
of the friction force (the starting friction force) which is e.ual to the product of the friction coefficient and 
the force normal to the sliding surface, During sliding the friction force is normally taken to be in the opposite 
direction to the movement, constant and equal to the starting friction force. Sometimes, the kinetic friction 
force is taken, according to Euler, to be smaller than the starting force, 


This paper deals with the problem of mathematical simulation of dry friction whose verbal description 
is given above, 


It is usual to base the mathematical simulation of dry friction on the formula expressing the friction 
force in terms of the sliding speed by means of Kronecker's formula 


T = —kPsignV, (1) 


where T is the friction force, k is the coefficient of friction, P is the normal force and V is the speed of 


sliding 
sign ¥ 1wihV>O, 


—1 with V < 0. 


Expression (1) does not represent the behavior of the friction force when the sliding contact surfaces are 
at rest. In a number of cases, however, when the starting friction force is equal to the kinetic, one expression 
(1) can be used both for the mathematical representation of the dry friction processes and for simulation, By 
using expression (1) it becomes possible formally to represent the forward rectilinear motion of the body by 
equation 


* Read at the seminar on the theory of automatic control at the Institute of Automation and Remote Control 
of the USSR, Academy of Sciences, 
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mi = F —kP sign x, (2) 


where m is the mass of the body, x is the coordinate determining the body's position, F is the external force 
acting along the friction surface, 


If |F| < kP, (0) =0, the solution k = 0 becomes a definite switching state of sliding [2] simulated by 
small high- frequency parasitic oscillations around the zero value of quantity x, This produces an effect similar 
to the effect of vibrational linearization; and the quantity x varies with a mean speed proportional to F, 


The additional verbal definition of function (1) as a function of two variables—the speed of sliding and 
the external force~is not used directly in simulating, but the discontinuous function of speed (1) is approxi- 
mated by a continuous one with a steeply sloping section, In this case it is possible to consider that dry friction 
at rest changes to large viscous friction, In order to simulate the excess of the starting friction force over the 
kinetic ong the approximating function of speed is varied so that its extremes come at the boundaries of the 
stated steeply-sloping section, In this case, as in the preceding one, the presence of a sloping section of the 
characteristic gives rise to parasitic deviations in the simulation process, 


The number of variations in simulating dry friction can be increased not only as the result of a more 
accurate interpretation of the dry friction law, by taking into consideration such circumstances as, for instance, 
the relation between the friction force and the duration of the state of rest contact, but also by examining 
various dynamic systems with friction for properties which could be eliminated. The latter method leads to 
relationships which cannot be obtained directly from simplifying equations of the type of [2] by neglecting 
certain terms, 


Let us now examine the forward rectilinear movement of a solid body with dry and viscous friction and 
a force of inertia, By using expression (1) let us write the motion equation in the form: 


mz = F —kP signa —lz —nz, (3) 
where 1 is the coefficient of viscous friction, n is the coefficient of the force of inertia, 


The representation of the static relation between x and F, which is attained at sufficiently slow varia- 
tions of F, by means of omitting the dynamic terms mx and 1x in equation (3), as it is done in[3, 4] leads 
to the relationship 


nz == F — kP sign z, (4) 


which does not represent the actual continuous relation between x and F, prevailing in this case, In the above 
mentioned papers formula (4) is corrected to a certain extent by additional verbal definitions, 


Let us also examine the forward rectilinear movement of a body with dry and viscous friction, but with 
negligible inertia, The relationship between the speed of the body and the external force for this case is 
shown in Fig, 1, a, If at the same time the excess of the starting friction force over the kinetic friction force 
is taken into account, the relationship shown in Fig, 1, b will be obtained, where the transition from the 
horizontal to the sloped section takes place at points F,; and —F, and the reverse action at points F, and —F, 
(5). 


Fig. 1 
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Let us now describe the results obtained by the author in simulating dry friction on model MN-8, 


The Case of the Starting Friction Force Being Equal to the 


Kinetic Friction Force 


As it has already been stated, the characteristic shown in Fig. 1, a represents the relations between the 
speed of the external force and the forward rectilinear movement of a solid body with dry and viscous friction 
and negligible inertia, This characteristic is produced with an opposite sign by the dead zone unit (DZU) of 
model MN-8, The mathematical representation of a more general case (taking into account the inertia of the 
body) can be obtained by adding to the external force the force of inertia, In simulating it is possible to use 
for this purpose the differentiating unit (DU) which differentiates and reverses the sign, 


~Q0909L 
a 
8 
~100N 
100 
b 


Fig, 2 


Moreover, if a movement with dry, but without viscous friction has to be represented, the viscous friction 
force with a reversed sign can be added to the external force, The circuit simulating the movement of the 
body in this case is represented in Fig, 2, a, where the block with the inscription VDU stands for the voltage 
dividing unit, The schematic of the DZU when used for reproducing the characteristic of Fig. 1,a* is given 
in Fig, 2,b, In Fig, 2,a the DZU is shown adjusted to remain insensitive up to 10 v and k the tangent of the 
slope angle of the external sections of its characteristic, is made equal to 1.1, The schegnatic of the block 
DU, which belongs to the auxiliary part of model MN-8, is shown in Fig, 2,c, The circuit contains small 
parasitic input resistances and a feedback capacitance, shown by a dotted line, When simulating friction, the 
input voltage to the differentiating unit was impressed through an additional capacitor of 0.1 uf, whose terminals 
were connected to the patch bay, With the additional capacity in circuit,the unit differentiated and multiplied 
by-0,0909 with parasitic time constants of 0.0024 sec and 0,0005 sec, The additional capacity was introduced 
to eliminate high- frequency oscillations, For the same purpose,differentiation and summation was carried out 
on high-frequency amplifiers model MN-8 with a band-pass of 1000 cps, The high-frequency amplifiers are 
marked on the circuit by the letter B, 


It is easy to check that the equation for the characteristic of Fig, 1,a has the form: 


* By means of the DZU of model MN-8 it is also possible to represent a characteristic with a limitation of 
coordinates, 
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By adding to the external force the force of inertia and subtracting the viscous friction force we obtain 
the equation: 


+ (F —c3X X + cg) (F— + ere + (6) 
2 . 


By substituting the coefficients with their numerical values according to Fig. 2,a we obtain 


— O12 + + — — 10), (7) 


It should be noted that by means of identical transformations of equation (7) (by multiplying its parts 
and the argument of the signum- function by a positive number and changing the variables) the absolute values 
of the coefficients can be chosen arbitrarily, Thus, it can be considered that the Fig. 2,a circuit, with the 
values of parameters shown in it, represents the general case of a forward rectilinear movement of a solid body 
with dry friction, Obviously, it is also possible directly to reset the parameter values of the simulating circuit. 
In some instances, however, this may produce high-frequency parasitic oscillations, Oscillations are caused, 
for instance, by the exclusion of the additional capacity of 0,1 yf, equivalent to increasing the gain of the 
differentiating amplifier (up to the parasitic high-frequency oscillations), 


A simulating circuit using differentiating units has a tendency to parasitic oscillation, It was not found 
possible, however, to solve equation (6) with respect to the higher derivative, and it was therefore necessary 
to use differentiating units when simulating dry friction in th’; case, For simulating in the differential form 
facilities are provided in the auxiliary part of model MN-8 for setting the initial conditions on the differentiating 
amplifiers, 


The results of simulating according to Fig, 2, a circuit are represented by oscillograms.* The oscillo- 
gram of Fig, 3,a represents a process which occurred in the simulating circuit and which corresponds to the 
following physical picture of movement, 


At first the external force F is equal to zerog next it assumes a value equal to half the starting friction 
force but the body remains at rest. Then force F suddenly jumps to a value one and a half times greater than 
the starting friction force, and sliding begins, The change in the sign of the external force leads to rapid 
braking and a change in the direction of sliding, When the external force becomes smaller than the kinetic 
friction force, the speed decreases and relative stability is restored, 


The only difference between the oscillograms shown in Figs, 3,b and 3, c and the one just described 
lies in the fact that in Fig. 3,b the VDU was not set to a gain of 0,909 but to 0,959 and in Fig, 3, c to 0.859. 
Also in the first case, the feedback coefficient of the DZU was greater than one,and in the second case smaller 
than one, which corresponds to introducing viscous friction with different signs, There is no distortion of the 
simulated processes in this case either, 


Above circuit can be used as part of another circuit for simulating more complex processes, Let us 
examine, as an example, the oscillations of an axis of a gyroscope freely suspended on a gimbal in a position 
nearly perpendicular between the suspension axis and the gyroscope rotor, with dry friction present in the 
bearings, 


Let us introduce the following notation: M, is the external moment about the axis of the gimbal internal 
suspension ring, y, is the angular velocity of the gyroscope rotor axis round the internal suspension-ring axis, 
I, the moment of inertia of the mass involved in movement round the internal suspension-ring axis, cy the 
moment of starting friction round the internal suspension-ring axis, Mg, y2, Iz, and cz are similar quantities for 
the movement round the external suspension-ring axis, and H is the kinetic moment of the gyroscope, 


* For the sake of a more convenient presentation the parts of the same curve at the points of discontinuity are 
joined by vertical lines, 
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Let us now write the equations of the moments about both axes of the gyroscope gimbal suspensions [6] 
by using the form of expression for dry friction, already obtained: 


= (My — 1191 — Hye + + — Fin + cays + 61) 


+ (My — Aye + — + ¢syi— ¢1) ‘ (8) 


+ (My + + — + sign(Ms Cas — 


Let us take the following numerical values for gyroscope parameters 1, = 5 g cm sec * 1,= 20 g cm sec’, 


Cy =c,=0.5 g cm, H=18,000 g cm sec, 


seals 


sec 


Fig. 3 


By giving their numerical values to the coefficients in equation (8) and using substitutions t = 300°! +, 
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z, = (20M, — 0,092, — 0,272, + 2, + 10) —sign 10) 


+ (20M, — 0,092, — 0,272» + z, — 10) 1+ sign (20Mi — 0.2722 + —10) , 


(9) 


= (20M, — 0,092, + 1,12; ++ + 10) 20Ms + 10) 


+ (20M, — 0,092, + 142, + 2, — 10) 1+ sign a + 1,42; + 22— 10) | 


Where the asterisk means differentiation with respect to rT. 


afl an pu 
V2 
2, 


Fig. 4 


The circuit of equations (9) set up on model MN-8 is given in Fig. 4, This circuit consists of two parts, 
each of which is equivalent to the circuit of Fig, 2,a and possesses all its properties. 


scale M, \M2 
(2112 127\190\@ 


Fig. 5 


The oscillogram of Fig. 5 shows the process of simulating the gyroscope movement about the internal 
axis of the gimbal suspension, due to a unit step action of M,. 
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It should be noted that in various actual gyroscopic instruments backlash in bearings and unbalanced 
suspension can considerably affect the gyroscope axis movement, 


The Case of the Starting Friction Force Being Greater Than The 


Kinetic Friction Force 


Let us now consider a more complicated case of dry friction when the rise of the starting friction force 
above the kinetic friction force is taken into account, We shall examine among the various simulating circuits 
not the simplest, which are attained with the help of so-called signature units of model MN-8, but those which 
are least subject to parasitic high-frequency oscillations, 


They can be examined in the same way as before if instead of the characteristic of Fig, 1,a that of 
Fig. 1,b is used, which can be represented formally as follows: 


+(F + + ys 
1 + sign (F + — ce) 1— sign (F + ys + ¢2) 


(10) 


’ ¥2 = — 


We shall ascribe to functions of the type y = sign (x + y) in their two- valued portion, both now and 
henceforth, the value equal to the preceding portion, which usually corresponds to the properties of the circuits 
embodying these functions, 


The first of the system (10) equations represents the relationship between k and the sum F + yy + yz 
similar to the relationship shown in Fig, 1,a, The second and third equations establish the relationship between 
the sum y, + yz and quantity F, a relationship shown in Fig. 6,a and realized by means of DZ units of model 
MN-8 according to the method shown in Fig, 6, b, 


I," % 


Fig. 6 


The simulating circuit of a body's forward rectilinear movement with dry friction, when the starting 
friction force exceeds the kinetic friction force, operating according to the characteristics of Fig. 1,b is shown 
in Fig, 6,c. 


The movement equations in this case will be 
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1 + sign (F — + + y2) 


= (F — Ce + + Yo) 2 
—si — + ee (11) 
+ sign t + — ¢2) 


1 — sign (F — + + yo + 
= — Cs 


By giving the coefficients their numerical values, which determine the parameters of Fig. 6,c circuit, 
we obtain 


x = 44,12) X 


2 


x + (1,1F—0,1X%+%+114+1,1y,+1,1y2) x 


4 —sign (1.17 —O,1z + 2+ 11 + 4.441 + 1,142) 


x 5 " 
1 +sign (1.17 + 2+ 1,1y, — 11) (12) 
{ — sign (1.1F —O 1x +2 + 1,142 + 11) 
¥2 = —9 2 


The oscillogram in Fig, 7 shows a simulating process in which the external force first changes from 
zero to 70% of the starting friction force, Next, due to a momentary rise of the external force above the 
starting friction force, sliding occurs which continues during the subsequent decrease of the external force 
until it reaches 70% of the value of the starting friction force (which is 40% higher than the kinetic friction 
force), The external force then changes the sign and as a result of braking a state of rest is restored, 


As an example of a more complex process whose 


geale 3 simulating circuit includes the one described above, let 
| atotaatity us examine again a gimbal gyroscope with a kinetic friction 


oe \ ie) ee force half the starting friction force and with the same re- 
maining parameters, The gyroscope movement equations 
in this case are easy to obtain from the gyroscope equation 

(8) by adding to M, and Mg quantities determined by re- 


£ pons lationships of the form of the second and third equations 
(11), In the simulation circuit shown in Fig. 8, the latter 


relationships were set-up by means of model MN-8 signature 
4 units, shown as rectangles extended in the vertical direction, 


A signature unit contains a polarized relay, which 
sec, shorts the middle output contact to the bottom one when 
Fig. 7 the unit input voltage is positive, and to the top one when 
‘ the voltage is negative. In order to decrease the parasitic 
processes when the relay is switched, the input voltage is impressed on the armature through an electronic 
trigger included in the unit, The input and the three output terminals of the signature unit are placed in the 
patch bay of the model, The input signals for the signature units are shaped by specially provided high gain 
amplifiers (for the purpose of decreasing the effect of the above mentioned parasitic processes due to relay 
switching), The feedback circuits of these amplifiers include diode limiters, not shown in the circuit, which 
cut-off the output voltage at 100 v thus keeping the amplifiers within the permissible limits of operation, 


The author intends to devote a separate paper to the use of signature units, 
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The oscillogram of Fig, 9, a shows the simulation process of gyroscope oscillations due to an external 
moment M,, about the internal suspension ring axis, a moment which is greater than the starting friction 
moment by 10%, The oscillogram of Fig, 9,b shows a process in which at first the external moments are 
equal to zero, then moment Mg of 0.3 g cm (the starting friction is 0.5 g cm) is applied about the external 
suspension ring axis. The system remains at rest, Next Mg assumes the value of 0,8 g cm and when move- 
ment begins Mg falls again to 0,3 g cm and the system reverts to a state of rest, With a repeated rise in My 
to 0.8 g cm the subsequent fall to 0,3 g cm occurs with a different type of gyroscope movement, Moreover, 
oscillations are produced at a lower frequency than in the case of the oscillogram of Fig. 9, a, 


The Case of the Starting Friction Force Depending on the Duration 


of the Contact At Rest 


Let us now assume as it was done in[7 and 8] that the excess of the starting friction force over the 
kinetic friction force is a function of the time the contact between the rubbing surfaces remains at rest and 
that it is expressed by formula 


T (2) = Too — (Too — To)", (13) 
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where To, is the friction force when the contact remains an infinitely long time at rest, Ty is the zero contact 
resting-time friction- force, equal to the sliding friction force, and r is the duration of the contact resting time. 


As before we shall base our reasoning on a relationship of the type of (5) whose parameter cy can be in- 
creased in relation to the contact resting time r. 
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When examining the forward rectilinear movement of a body with dry and viscous friction and insignificant 
inertia, let us use the approximate relationships 


} 


+ sign (F—ce— 


Yi = ys c2—y1) 1 + sign + y1) (14) 


ign (F — — sign (F 
a + ya {es + 2 v1) sign 


1+ sign (F + ce + 1 — sign (F — 
2 2 


These equations represent the process with increasing accuracy as cg increases, It is easy to see that the 
right hand side of the third equation (14) cannot be negative, With a positive initial value of yz its subsequent 
values will also be positive, Quantity yy which, according to the second equation, is equal either to zero or to 
yz Will not be negative either, Thus the value of the friction force, determined by the first equation (14), has 
a minimum equal to cy and can rise on account of the y, term, 


At first let |F| < cp+ yy expressing a stationary contact, Equation (14) will then take the form: 


= (F — C3 — sign C2 4. (F Cs sign ( st + Yi) : 
Yi = Yo» Yo + = C5. 
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Equation (15) represents the rise in the starting friction according to relationship (13), Now as a result 
of external action let |F | > cg+ yy. In this case equations (14) become equivalent to equations 


= (F — eg — y,) + yy) 
(16) 
yi = 9, Yo + (C3 + C4) Yo = 9, 


Equations (16) represent sliding with friction cz, Moreover, the larger the sum cs + c,, the quicker will 
the intermediate variable y, tend to zero, If it happens further that,as a result of an external action, | F | 
becomes smaller than cg, the process will then be represented by Equation (15). 


D B A Let us note that from equations (14) it is possible to pass to less 
Ee Vv accurate but at the same time simpler equations, for instance, by 
c substituting y, for y, in the first and third equations and by omitting 
the second equation, In addition, the right-hand side of the third 
Fig. 10 Ane can be substituted by cg providing the value of cy, is sufficiently 
rge. 


Figure 10 shows a device in which cylinder A, moving with a constant velocity V, draws by means of 
spring B a body D over surface C, Under the influence of the rise of the starting friction force above the 
kinetic friction force there arises an intermittent movement of the body, 


In simulating above case of dry friction,equation (14) was used, somewhat simplified in the manner 
indicated above at the expense of the right-hand side of the third equation, By adding to the external force 
the force of inertia and subtracting the viscous friction force, and by assuming the external force to be de- 
termined according to Fig, 10 we obtain equation 


= (F — ct + — — ys) + 4 


(17) 


_ — sign (F — eet + cx — cg ~ yx) 1 + sign (F — cot + + + 1) 
Yi = Ye 2 2 


2 
es + y1) } 
F=V—z. 


The simulation circuit corresponding to equation (17) with cy =1, cg= 10, cg =0,2, cg = 500, cg = 3, 
Cg = 0.0909 is shown in Fig, 11, In it,signature units containing two polarized relays, operated by the same 


input signal, were used, Both in the circuit and the patch bay,the three top output contacts are switched by 
one, and the three bottom ones by the other relay. 


The oscillogram of Fig, 12 represents the process of simulating intermittent motion arising, after a long 
state of rest, as the result of a unit step action V, It clearly shows the peculiarity of this process consisting in 


the first speed pulse being larger than the rest, This peculiarity was first explained by A, Iu, Ishlinskii and 
I, V. Kragel'skii [7 and 8), 


It should be noted that in practice, when simulating according to Fig, 11, the drop in the intermediate 
value of y, at the beginning of the simulation of sliding can develop slower than it should according to equa- 
tion (17) since the voltage at the integrator input, shaped by the summing amplifier, is limited by the diodes 
to the scale of the model, i.e,, to 100 v. In this case it did not influence the results of simulation since the 
duration of the speed pulses in the process under consideration was sufficiently large. The parasitic discharge 
time of the integrating condenser can be reduced to a minimum by appropriately shorting its plates through 
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a signature unit, This possibility, however, is considered to be beyond the scope of the present investigations 
since the patch bay of model MN-8 does not provide for such an operation, 


Fig. 11 
The above results, obviously, do not exhaust the problem of simulating dry friction, Its examination 
can be extended, For instance, it is useful to examine the simulation of the forward movement of a body over 
a plane, or the movement of a shaft tightly seated in a cylindrical bearing with two degrees of freedom, and 
similar but more complex cases of nonlinear processes involving dry friction, 


_ scale 


Fig, 12 


In the first of the above cases the motion equations can be written as follows: 
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(Fx —¢,X + +(Fy—ea 
Ye = yi Signy, 
1 + sign (ys — ¢1) 


Ys = (Y2 — 2 ’ 
F,— Cat + 
Yo 3 
Ey — cay + cay 
Coy = Ys > 


where cy > Cs => 0, Fy and Fy are the components of the ex- 
ternal force along the axes of Cartesian coordinates, y;, yx, 
and ys are intermediate variables,and x and y are the 
coordinates of the moving body, 


The starting friction is taken to be equal to the kinetic 
friction and is represented by c,. 


The movement of a shaft, tightly seated in cylindrical 
bearings, with two degrees of freedom, that of rotation and 
sliding along the axis, when the difference between the 
starting friction force and the kinetic friction force is taken 
into account, can be expressed by the system 


— + 4 t= ey + (19) 


cy C2 


Yo = Signy, 


, 1 
ys = sign(ys +: ¥s— +*), 
1+ si 
Ys = (Y2 — Cs) ’ 
Y2 + Cg 5 
M — evy + esy 
Coy = Ya i—si 
+ Cg sven 


where 1> cg > cg>0, F is the external force, M is the external moment, y;, ya, ys, and yg are intermediate 
variables, x is the longitudinal displacement of the shaft, y is the shaft’s turning angle, cy is the starting friction 
force, cg is the starting friction moment, cz is the fractional expression for the kinetic friction force in terms of 
the starting friction force, 


The general case of a heavy homogeneous disc sliding along a plane with dry friction (with the starting 
friction force equal to the kinetic friction force) can be represented by the system 


F,—mz+ 


F,—my j 
y = V (4, 29) ¥ (a0) 


21 
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= (24, 23), 
= V (Fx — mz + + (Fy — my + 
2, = M — la + cya, 


where functions V and w of the two variables z; and zz have the form shown in Fig. 13 and the square root is 
used in its arithmetic meaning, Also where x and y are Cartesian coordinates of the center of the disc, a 
is the turning angle of the disc, Fy and Fy are the components of the principal vector of the forces acting on 
the disc, M is the principal moment acting on the disc with respect to its axis of rotation, m is the mass of 
the disc, I is the polar moment of inertia of the disc, 


When functions V and w were plotted, in addition to dry friction, viscous friction at the points of contact 
was also introduced and then the inverse functions — z,(V; w) and + 2z4(V; w) were calculated, these functions 
having the meaning of the principal moment and the modulus of the principal vector of the dry and viscous 
friction forces with respect to the disc axis, By selecting suitable values for coefficients cy and cy the effect 
of viscous friction on the final result. can be eliminated in the same way as it was done previously, 


SUMMARY 


1, Differential equations for sliding bodies with dry friction under various conditions have been derived, 
The equations thus derived are distinguished by the absence of the usual verbal additional definitions, the 
absence of sliding solutions, and in addition they are not solved with respect to higher order derivatives, 


2, Simulation of sliding with dry friction was carried out according to equations which were not solved 
with respect to higher order derivatives, by means of differentiating amplifiers provided for the purpose in 


3, When simulating by means of differentiating amplifiers, the initial conditions, which determine the 
particular solutions of the equations, are introduced by charging the differentiating condensers in a manner 

similar to the one used in simulating by means of integrating amplifiers, Model MN-8 is the first of its kind 
using this arrangement, 
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model MN-8, 


CHOICE OF A POWER UNIT FOR AN OPTIMUM AUTOMATIC CONTROL SYSTEM 


L. N, Fitsner 


(Moscow) 


A method of choosing a power unit for an optimum automatic control system is 
recommended, 


To illustrate the recommended method a servosystem whose power unit consists 
of a dc motor is given as an example, 


1, Formulating the Problem 


A number of articles [1-10] have appeared in recent years dealing with optimum automatic control 
systems (ACS), 


A common posing of the problem is characteristic for all these papers, They take one part of the ACS 
for granted (usually the power unit) and examine the design of a control system with a minimum response time, 


The successful solution of this problem now provides opportunities for a more correct and valid approach 
to the choice of a power unit as well, 


In fact for an optimum ACS with given limitations and a fixed power unit and loading there exists a 
definite relation between the input signal law of change and the response time, 


In this connection it is expedient to know this relation in advance for the most common types of power 
units and loading in systems working with some typical, for instance, step-by-step disturbance, 


If such information is obtained and presented in the form of a family of curves, it will provide the de- 
signer of AC systems with answers to a number of important questions, 


With their help it will be possible in the first place to determine the response time of an ACS with any 
given power unit, 


If the response time, obtained from the curves, should exceed the permissible value, one can immediately 
conclude that the preliminary choice of the power unit was not suitable, since even in a system with an optimum 
control it cannot provide the required speed of operation, 


It can happen that with a given load the system will not be able to ensure the required speed of operation 
even with the best chosen power unit, 


The curves will provide an answer to this most important question as well, 


It will also be possible to tell from these graphs how near to the optimum is an ACS with a nonoptimum 
control unit, in order to be able to decide whether it is expedient to change the control unit for the purpose of 
improving its time characteristics, 


The technique of plotting such curves is shown by an example given below, 
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The design of a servosystem with optimum operation is also examined, its time characteristics, obtained 
experimentally, are compared with those obtained from the graph, 


A simple relay type servosystem powered by a dc motor with separate excitation is given as an example, 
The motor operates a platform with a constant moment of inertia. 


The restricted value of the system is the voltage at the motor windings, 


2. Motion Equation of an Optimum ACS 


Figure 1 gives the schematic of a servosystem, 
Motor M, with a separate excitation, operates platform Pl, 


The control unit Ctl calculates the best control law for the motor and changes the veltage polarity at 
the armature windings by means of relay P contacts, 


16 
\ 
+o 
Pot 
Pl 
7 
re 
Fig. 1 


The platform movement is proportional to voltage Ujp, which is impressed together with its derivative 
dU; ,,/dt on the control unit input, The voltage from potentiometer Pot slide (U4,,) and tachometer generator 
TG (dUg,,/dt), proportional to the platform deviation angle and the platform speed, respectively, are also im- 
pressed on the input, 


A second relay P, can be provided at the output of the control unit for connecting the motor armature 
winding, after balancing, to the linear amplifier, or if required for disconnecting the winding, 


The transient process in the motor can be represented by the following equations: 


U = iR4-L + co, (1) 
My, = Col, (2) 

dw 
My =J—-, (3) 


where U is the restricted in value direct voltage at the armature winding, i is the armature current, R is 

the armature resistance, c; is the back emf coefficient, w is the angular velocity, cz is a constant coefficient, 
J is the moment of inertia of the platform and motor, referred to the motor shaft, L is the armature inductance, 
and Mj, is the motor turning moment. 


In writing these equations certain factors were not taken into consideration, such as the friction force 
in the gear train, armature reaction, eddy currents and hysteresis in the motor, 


The armature self-induction coefficient was assumed constant, 
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Experimental data, given below, shows that for the case under consideration the omission of above 
quantities did not make any appreciable difference to the duration of the transient process, 


Using (1) - (3) let us write the velocity differential equation in the form: 


dw» Rdw = 
1 da 1 | U 
8 


where T, = is the electrical time constant, = 
R CyCe 


The velocity equation for dynamic braking by reversing the armature will differ from (4) only by the 
sign of the right-hand side of the equation: 


1 dw 1 ae, 


Let us henceforth denote respectively by x;; and wy, the motor turning angle and its velocity at the 
beginning of the movement before braking starts; and by xy and wy the corresponding quantities for the 
second, final interval during dynamic braking, 


Let us first examine the transient process when a step voltage is connected to the system input at instant 
t=0. 


Let us assume that at the moment the signal was applied the system was in the state of rest, i,e,, w ;=0 


and Se. =0 att=0, The solution of equation (4") under these conditions will take the form 


dings pal prU U (6) 
(Pi--P2)er (Para) 
where 
1 i i (7) 
27 e 4T% le Tom 
i (8) 


The motor turning angle at the starting interval will be 


tn 
Z_ = ou de. (9) 


0 
Where ty; denotes the starting interval duration, By using equation (6), we obtain 


Pil Pi — Po(Pi — Pa)er o, Pi( Pi — Paden Pol Pi-— Paer (10) 


= 


We find the turning angle and motor speed in the second interval by integrating (5), 


Let us denote the motor shaft angle and motor speed deviations from the required values at instant ty, 
when the voltage polarity at the armature is reversed, by 2, and @,. respectively, We shall then obtain 


* Below we shall only examine systems for which py and pz are real quantities (T,,, = 4T¢). 
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for the starting instant of the second interval when t =0, 


U ) Pr U U (11) 
= 
‘x (12) 
= onde, 


0 
Where ty denotes the duration of the second interval, By integrating (12), we obtain: 


P2(P2 — pi) 


The obtained solutions of the equations will be required below for deriving the system time characteristics, 


3. Control Unit Block Schematic 


Let us now examine the arrangement of an optimum system control unit on whose input is impressed, as 
distinct from the previously examined case, a signal of frequency wy. 


In the first place let us find the relationship between the velocity and position balances, at which the 
control unit must reverse by means of relay Py the voltage polarity at the armature, 


In the second interval when dynamic braking is applied the system will deviate by angle Xx equal to 


= + Opty + (14) 


The second right-hand side term of this expression accounts for the change of the input signal during 
the second interval t,, The rate of change of the input signal ws over a small interval of time ty is assumed 
to be constant, 


Equation (14) third right-hand side term accounts for relay P (Fig. 1) operation time top, It is often 
impossible to ignore this quantity since it is of the same order as the transient state duration, 


Considering that when t=ty, w= ws, we obtain 


U\ Pi ( OV 
With the help of (13) and (14), we obtain 
U U 
Ta Mate Ont = — — Pi) Pa (Pa Pi) 
U 
U Pi(w,+ (16) 


In this case it is impossible to find the required relationship between x and wy by the method 
suggested in [1] since velocity ws at the end of the interval is not zero, Hence, for obtaining this relationship 

from (15) and (16) let us utilize a computer circuit shown in Fig. 2, 


This circuit controls relays P and P, (Fig. 1) according to the input signals which are proportional to the 
required and actual values of the velocity and deviation angle, The circuit consists of four operational ampli- 
fiers ¥,; - Y4[11], used for amplifying, adding and measuring the polarity of the input voltages, of two multiplying 


* This equality is valid because at the instant the relay contacts are switched the current stops flowing through 
the armature, It is shown in [10] that with the adopted method of switching the transient state duration is de- 
creased, For the same reason the optimum process under consideration in this system is represented by two in- 
tegrals only,although the system contains a third order restriction, 
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links LML and 2ML [12] which multiply the input voltages, and of two nonlinear transducers INT and 2NT 
[13] whose output voltages are nonlinear functions of their input voltages,* The values of the voltages in 
various parts of the circuit are denoted by the same letters as the variables in equations (15) and (16), Thus, 


a denotes the sum = 
| 


eat 
ent 


Fig. 2 


At the output of amplifier Y, there appears a voltage proportional to the minimum time required for 
making w = wg, The term for correcting the error due to the relay operation time wtpp appears at the output 
of amplifier Ys as the result of multiplying velocity w by the amplifier Yg fixed gain, equal to the relay 


operating time top: In order to make this circuit control the motor two relays P and P; are connected to the 
circuit output, 


In certain cases the electromechanical time constant Te,, ‘s considerably larger than the electrical one 
Te, and the latter can be neglected, ** The computer control unit circuit then becomes considerably simpler, 


At Te =0 instead of (15) and (16), we obtain 


f Ta. = U 
tx 


With the help of (17) and (18) we obtain 


U U \ 
U (19) 
— sign In ( Wy | +- 
Where sign ws and sign wy, show that with the change of the signs of ws and wy, the signs of 


p In (| @, | + =) In + =) change accordingly. 
1 
e The circuit of the optimum system control unit for the case in question is shown in Fig, 3, It includes 
.- five operational amplifiers (Y, - Y,), two nonlinear transducers INT and 2NT one multiplying link ML and two 
plying 


* The method of obtaining circuits for solving transcendental equations is examined in [14], 
** The possibility of neglecting T, can be judged from a graph given below, 


}—-{2ML| 
igh 
n- 
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relays P and P, for switching the armature, 


At the output of the first nonlinear transducer there appears a voltage proportional to sign @, In ( 3 | + —), 
1 


and at the output of the second one a voltage proportional to sign w, In ( ny | + — Such relations can be 
1 


obtained fairly simply by means of carborundum resistors (tyrite and vilite),* 


The first and last terms of the figure brackets of expression (19) are represented in the form of a product 


> Saeed which is formed in the computer by one of the resistors at the input of Ys. 
em 


The cperational amplifiers Y; = Y, are shown in Fig. 3 only for the purpose of making it simpler to follow 
the calculation carried out according to (19), 


This circuit can be simplified considerably and 


4, need not contain above amplifiers; their functions 
Be} can be taken over by the nonlinear links LNT, 2NT, 
@c = on and ML. 
a 4 
- + Finally, if it is assumed that the system operates 
Fig. 4 3 on a step-by-step unbalance only, w3=0 and the values 
of the velocity and angle unbalance will be related by 
the expression 
c,@ 


or taking into consideration the relay operating time 


U 
tu = Em 


The control unit circuit for this case is shown in Fig, 4, It includes amplifier Y,, tyrite Ry, resistor R2, 
connected in series with the tyrite and resistors Ry and Rg, Resistor R, is connected in series with the tyrite 
in order to make its nonlinear characteristic close to the required one, Resistor Ry is included in order to com- 
pensate the effect produced by the relays P and P; operating time top on the control signal, 


These relays, in a manner similar to the preceding circuits, are meant to control the servo motor. 


* Carborundum resistors are also used for an exact and simple representation of the relation y= x |x|, required 
for optimum systems examined in [1], 


re 
| Fig, 3 
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4. Graph for Choosing a Power Unit 


The design of servosystems is usually started by choosing a power unit, For a correct choice of a motor 
(or any other type of power unit) it is necessary to know the relation between the load, the motor parameters 
and the response to at least one typical, for instance, step-by-step input signal, 


, oo It was stated above that such relationships could be ob- 
7 ee “nu ; tained for certain types of optimum servo systems and represented, 
for instance, in the form of a family of curves, 


a 


mx Such curves for the type of servo systems under considera- 
} + tion are shown in Fig. 5, They have been plotted on the basis 

' of equations (10) and (13) for the two intervals of the system 
y movement considered, 


XT 


Baeeee ae For a more convenient graphic representation the above 
3 7 relations are given in Fig, 5 in dimensionless units, Along the 
axis was taken the dimensionless response time 


Fig. t 
(22) 


and along the Y-axis the dimensionless quantity 
re 
(23) 
em 
Where t_ is the system response time to a given unbalance and x is the value of the unbalance, 


As a parameter for the family of curves serves the quantity 


(24) 


The curves for n = 5 and n = are close to each other, Hence it is possible to neglect the electrical 
time constant in this case (roots p; and p, are real), Moreover, as it was pointed out before, the control unit 
circuit is thus considerably simplified (Figs, 3 and 4), 


It is possible by means of these curves, on the basis of the given platform moment of inertia and the per- 
missible response time.to choose a required type of motor or to arrive at the conclusion that with the given 
servo system it is impossible to meet the given technical requirements, 


If the type of motor is given,it is possible to determine the response time for the required platform 
moment of inertia. 


By means of the same curves it is possible to determine how close to an optimum is any servosystem. 
Such comparisons must, of cours., be made with similar loading, As the result of these comparisons ‘t is 
possible to decide whether it is expedient to change the control unit of any given servosystem with a view 
to improving its speed of operation, 


5. Experimental Investigation 


A most simple optimum servo system, whose circuit is shown in Fig, 4, was investigated experimentally. 


The investigation was carried out with a view of checking its efficiency and for the purpose of evaluating 
the error due to the idealization of the system in its mathematical treatment,* 


The model of an optimum servosystem included a dc motor type SL- 267, which turned, by means of a 
reduction gear of ratio ip = 120,a platform with a moment of inertia J = 0,05 kg sec’, 


* V. P, Grekova took part in the experimental investigation of the model of the system. 
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A voltage proportional to the velocity of the platform was obtained from tachometer generator whose 
shaft was connected by means of reduction gear to the motor shaft, 


A voltage proportional to the turning angle of the platform was 
a obtained fromm the slide of a wire wound potentiometer fixed to the 
platform shaft, 


mentally at T,,,, = 0.1 sec by changing the moment of inertia of the 
platform, A voltage of 110 v was connected to the armature winding, 
The nonlinear ratio U/c,, required for calculating the system re- 
sistance, and equal to the stable state yee | of the platform, was 
also determined experimentally at 3,76 sec “. 


as 
20 The electromechanical time constant was determined experi- 


The voltage at the potentiometer slide (Uy) was proportional 
to the turning angle of the platform 


=> kyx, (25) 
where ky = 2 v/degree, 


26} The voltage Ug at the tachometer generator was proportional 
to the velocity of the platform: 


(26) 


Urg= ’ 
where k, = 0,083 v sec/degree, 


The switching of the armature was made by the contacts of relay Py whose winding was connected to the 
output of the operational amplifier Y, (Fig. 4), The relay operating time was equal to top = 20 sec, 


Figure 6 shows the required relation between the velocity and position unbalance at the instant of the 
armature voltage switching (curve 1), Voltage U proportional to the velocity was plotted along the X-axis 
and current i proportional to the right-hand side of formula (20) was plotted along the Y-axis, The same 
graph shows the volt-ampere characteristic of a tyrite (curve 2), The tyrite was in the shape of a cylinder 
10 mm high and 50 mm in diameter, 


In order to obtain a better agreement between the natural volt- ampere characteristic of a tyrite and 
the required curve, resistance Rg (Fig. 4) was connected in series with the tyrite, The points of the volt- 
ampere characteristic of such a combination are shown in Fig. 6, Practically all of them coincide with the 
required curve, The value of Rs can be chosen on the basis of the following considerations, 


Let the switching occur at the maximum velocity wy = Wmax = = » Which corresponds according to 
(20) to the position unbalance of 


(27) 


Xn = Tem Wmax In 2) = 0.31 Tem“max: 


With selected values of Tey, and Wmax = 0,116 radians. According to (25) Ug =13.3v and 
according to (26) =17.9 v. 


At the instant of contact switching the voltage at the amplifier input (at point a, Fig, 4) is equal to 
zero, Hence, 


(28) 


U 
According to Fig, 6 = 24,5+10°° amp at wax. Thus, Rs =544 kohms, It is possible to 


find the value of Ry which compensates for the relay operating time from (20), (25), and (26): 


4 
Fig. 6 
|| 
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With the chosen circuit parameters R, = 1,1 megohm, 


Figure 7 shows oscillograms of the response time of the optimum servosystem under consideration, 
working with step unbalances of 20, 30, and 40° (oscillograms a , b, and c, respectively). 


The oscillograms show the values of the posi- 
tion (x) and velocity (w) unbalances and the voltage 
4 at the amplifier Y (U) output, 


1 = At the end of the transition process the position 
‘s es as and velocity unbalances became equal to zero, The 
a1 switching of the relay contacts occurred later than 
Ot the jump of the amplifier Y output voltage by the 

w amount Of top. 


| \ When tested the system did not include the 


linear amplifier shown in the circuit of Fig, 1. 
Relay P; (Fig. 1), which disconnects the motor arma- 
ture at small unbalances,was also missing, That is 
dM sec > why at the end of the transition process oscillations 


T were observed, 
° ae Let us determine from Fig, 5.the response time 
— 


b of the system and compare it with that obtained from 
w the oscillograms. 
According to (23) = = 0.93 atx= 
Q/ sec = 20°, This value of y corresponds to r = 2,06 and 


t® = rTe = 0.205 sec, 


8 


The relative error in determining the response 


' Beh sabe a time by the graph of Fig. 5 as compared with the 
z 
c experiment is equalto 6 = 10%, 


w Above response time was obtained from the 
wl A oscillogram of Fig, 7, a, For oscillograms in Fig. 7, 
oT b and c similar errors were found to be equal to 
Fig. 7 8.9% and 8%, respectively, 


Thus, for the three cases examined the actual response time did not differ from the calculated one by 
more than 10%, 


If the time constant T,,, should change in the course of the operation of the system (for instance with 
time or with the deviation angle), the optimum law of control can be maintained by an automatic variation of 
resistances R, and Rg in the control unit, 


SUMMARY 


It is expedient to plot curves showing the relation between the speed of operation of an optimum AC 
system and its power unit parameters for the most common types of power units and loading, 


It is possible by means of these curves to find the response time of an optimum ACS which includes any 
type of power unit, thus, determining the usefulness of the latter, 
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For the type of optimum systems investigated the response time obtained experimentally differs by not 
more than 10% from the one found from the curves, 
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CALCULATION OF TIME CHARACTERISTICS OF 
PNEUMATIC FLOW CHAMBERS 


Vv. N. Dmitriev and V, I, Chernyshev 
(Moscow) 
Nonlinear differential equations of pressure variations in a pneumatic flow 
chamber are integrated with the assumption that any of the input quantities are 


subject to step variations,* Examples of time characteristics calculations are given 
and these results are compared with characteristics obtained experimentally, 


1, Equations for Pneumatic Flow Chamber Time Characteristics 


The adopted notations are: Q* is the weight of air in the interthrottle chamber, V is the volume of the 
interthrottle chamber, y * is the specific gravity of air in the interthrottle chamber, R is the gas constant, 
T is the absolute temperature of air, G*, is the weight rate of flow of air through the i-th throttle, P, is the 
feeding pressure of the pneumatic flow chamber, P*%, is air pressure in the interthrottle chamber, Pz is the air 
pressure passed the second throttle, f; = 4F, the effective port area of the i-th throttle, yj; is the rate of 
flow coefficient of the i-th throttle, Fj is the port area of the i-th throttle, i is the number of the throttle, 


andt is time, 
Let us now deduce a differential equation for pressure Py* variations 
} t with respect to time in the pneumatic flow chamber (Fig, 1), Let us 
| assume that Ty =T,=T,=T. 
c 
4 ES iis The weight of air in the interthrottle chamber will be 
Fig, 1 


Differentiating equation (1) with respect to time and using the gas equation we obtain 


(2) 


In equation (2) quantity dQ*/dt represents the weight rate of flow of air into the interthrottle chamber 
(or out of it) in a dynamic state, i,e,, dQ*/dt =Ge —Gz, 


Taking this into consideration and rewriting equation (2), we obtain 


aP; 


(3) 
RF — Ge 


The weight rates of flow Gy* and G,* can be calculated from formulas [2]: 


* Differential equations of pressure variations in pneumatic device flow chambers with linear approximation 
are considered in [1], 
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a) for a pre-critical state of flow 


P 
Gi=hy/ Rr Pi(Pia—Pi), 0,5: 


b) for an above-critical state of flow 


P 
Gi = ur < 0,5. (4b) 


It should be noted that there can be four possible combinations of flow through the first and second 
throttles: “a-—p", a" and “a—a® (pre-critical flow through the first and second throttles, 
above-critical flow through the first and pre- critical through the second, etc,), 


By substituting in (3) the corresponding expressions in (4a) and (4b) it is possible to obtain four differential 
equations related to the four combinations of flow. The expressions thus obtained, subject to the condition 
that any of the input quantities (Py, f;, f2, and P,) changes in a step-by-step manner, become equations of 
the first order with separable variables and can be integrated, In integrating let us assume that the rate of 
flow coefficients jy are constant during the transient state, As the result of integration, we obtain: 


for combination “p — p* 


A 2r,;—1 (2m. — 1) +r: 2r + | 
(it—ri) 29 281 Ba V ri (4 — ri) + B*)(r,—r) 


BsBa 
In 


V 


where c is the integration constant which is found from the initial conditions, 


P. P . 4 


for combination “a— p" 


ry = P,/P3; 


where 


for combination a" 


2 VY r 
t= [2k are tg (7) 
and for combination “a — a” 
t=Aln|k—r,| +c. (8) 


* Arc tg=tan In =log.— Publisher's note, 
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| 
* 


In the case when f, or f, become zero, expressions (5)- (8) degenerate respectively into equations re- 


presenting the emptying and filling of the dead-end chamber, 


root 
= 
T 
a6 “SSN as 
NY 
a5 
a \ 
a2 
a 
as 1 45 2 
Fig. 2 


Case when f; =0: 
a) pre-critical flow 


b) above-critical flow 


Case when f,=0: 
2) pre-critical flow 


t= Aarctg 


ri + 
where 
V 2 
A= 
b) above-critical flow 
t => Ayr, +- Cc. 


(9) 


(10) 


(11) 


(12) 


It is interesting to note that the transient processes in the interthrottle chambers do not depend on the 
absolute values of Pg, Py* , and Pz, but on their ratios r, rz, andr. For calculating the transient process in the 
interthrottle chamber it is necessary to know beforehand through what flow combinations the process will pass 
(and it can, as we shall see below, extend to three flow combinations through the first and second throttles), 
and to what state of flow the initial static state corresponds and what are its parameters, It is impossible to 
solve these problems by means of static equations for a pneumatic flow chamber, since if Py* is unknown, the 


combination of flows through the first and second throttles is also unknown, and it is impossi vie to decide which 


of the static equations should be used (corresponding to “p — p", “p — a", “a —p", or “a ~ a® combinations), 


~-—i +e; 
re 
1 
t= —Aln—+e. 
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In addition when calculating transient processes in the inter- 
throttle chamber of a pneumatic device (Fig, 1) it is necessary to 
know the rate of flow coefficients 1; through the throttle devices, 
A pneumatic flow chamber in an auto- pneumatic devices usually 
takes the form of a nozzle-vane type relay (Fig. 3), The relation- 
ship between the rate of flow coefficient and the throttle para- 
meters for such a relay are usually given in a form of a graph, 

A graph representing this relationship for a nozzle-vane type 
throttle is shown in Fig, 4, A description of two such graphs is 
given below, 


NN 


ye ke 


a02 006 AN 020 0% o2eh,mm 


2, Auxiliary Graphs 


Figure 2 shows a graph by means of which it is possible to determine the initial static state parameters, 
the combination of flows at the initial static state, the combination of flows through the first and second throttles 
during the transient process and the combination of flows at the new stable state value of the variable parameter, 
The graph is split-up into four zones, corresponding to the four possible combinations of flow through the first 
and second throttles, 


The equations shown in Fig, 2 are obtained from the condition of equality of static weight flows [ex- 
pressions (4a) and (4b)} through the first and second throttles for various flow combinations, These equations 
have the form: 


for the region of combinations “p — p* 


(2) = ry (13) 


for the region of combinations “p — a" 


* A similar graph for Saint-Venant-Wentzel formulas is given in [3], 
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: This problem is solved by means of curves in Fig, 2 plotted from static equations for a pneumatic flow chamber,* 
i 3 
Fig. 3 
| 
Fig, 4 
| 
for the region of combinations “a — 
(ty = 4r? 


eter. 


for the region of combinations “a — a” 


An equation for the boundary between regions of flow combinations “p — p" and “a — p" can be found 
if in equation (13) r is substituted by mrp and it is assumed that ry =0.5, After this transformation we 


obtain >) =1_(1 — ry). 


The equation for the boundary between the flow regions “a— p", “a— a" and “p—p", “p—a™ is 


r2=9,5 and between regions “a — a" and “p—a™ it is = 0.5, 
2 


It will be seen from Fig, 2 that in order to find pressure P, from the graph the ratio of effective areas 


7 and the ratio of pressures r= _ must be known, For a given value of r and $- , it is possible 
to find from the graph the value of ry = P»/P;, from which with the knowledge of P, it is possible to find P,, 
This also determines the flow combination which corresponds to the initial static state, It is possible to de- 
termine from the graph the flow combinations covered by the transitional state, Thus, with r= P,/P) = const 
and a step-by-step changing port area one should follow the given r curve from the point of the initial static 
state to the point of the new stable.state value, 


Thus, the regions covered by the transient characteristic can be found and,hence,the equations which 
should be used for its calculation, The points of intersection of curve r= const with the region boundary lines 
will provide the boundary values for quantities rz and ry. 


In a similar manner it is possible to find the air flow combinations through the throttling elements with 
a step-by-step variation of Py and P, at f/f, =const, with the one difference that in this case the various 
flow-state regions should be traversed vertically, 


Figure 4 shows the relationship between the rate of flow coefficient yj, the distance h from the nozzle 
to the vane and the pressure Pj. in front of the nozzle of the nozzle- vane type throttle, The rate of flow 
coefficients 4; have been determined experimentally [4]* and include a range of nozzle diameters from d, = 
=0.3 mm to d,;=3 mm at T = 239° K, P2=1 A,*and D/d,=1,6 to 6.5, 


3. Examples 


In order to illustrate the technique of calculating the time characteristics of pneumatic relays of the 
type under consideration by means of the obtained formulas and graphs we shall give four examples, An 
actual pneumatic relay of the nozzle-vane type with the dimensions shown in Fig, 3 was taken for an example, 
Time characteristics, with a step-by-step variation of the second throttle device port area and the feeding 
pressure Py, were obtained experimentally for this relay, thanks to which it was possible to compare the cal- 
culated and experimental time characteristics, 


Example 1, The given data was: Py) =2,5 As, P2=1 A,,* T = 293° K (Fig, 3), The distance between 
the nozzle and the vane changes in a step of hy = 0,037 mm to hgrapje = 0,083 mm, It is required to derive 
a time characteristic of the pneumatic relay P, = f(t), 


Before attempting to derive the transitional process it is necessary to determine the pressure Py in the 
interthrottle chamber at the initial static state and the pressure Pj stable at the new stable state condition, In 
connection with the fact that the rate of flow coefficient, for the second throttle device is not constant and 
depends with a constant opening h from the pressure in front of and behind the nozzle, the values of the above 
pressures should be obtained by the method of consecutive approximations, 


* The dotted line has been obtained by extrapolation, 
** Here and hereinafter: A, = absolute pressure, 


1091 


he 4(re—r)' 
fe 
ttles | 
3 


1092 


Let us find the value of Py. Taking into consideration that at the initial static state hy = 0,037 mm and 
the pressure in the interthrottle chamber lies between the limits of 2,5 A, = Py = 1 Ag,let us take as the first 
approximation 3 = 0,4 (Fig. 4), The constant throttle consists of a hole in thin wall, hence it can be assumed 
that 1 = const = 0,8 [4], 


Let us find the ratio 


2 


and by means of the graph in Fig, 2, considering that r=0.4, find r2= 0.47 and hence Pyp = 2,12 Ag, 


From the graph Fig, 4 we find a more accurate value for ue =0,48, Then (f;/ f ~ ss 1 and hence 
Pio =2 Age 


Let us find the stable value of Pystapje. As 4 first approximation for hstah = 0.083 mm we find from 
graph Fig, 4 that ut =0.8 and determine (f;/ f,)! =0.218, By means of graph Fig. 2 we find rz = 0.9 and 
Pistab = 1-08 As. Using the graph Fig, 3 again we obtain a more accurate value for u # = 0,755 and cal- 
culate (f,/ f,) = 0.29. 


From graph Fig, 2 we determine ry = 0.877 and hence Py = 1.14 Age 


The value of k is equal to 0,29, i,e., to the ratio f;/ f; in the newly established stable state, Let us 
also calculate f, for the new state, remembering that 2=0.755; f2= 2 = 0,004 cm’, 


Now we have all the data required for calculations, It only remains to determine the flow combinations 
of the first and second throttles covered by the transient state, For this purpose let us use the graph of Fig, 2, 
Marking off on curve r =0,4 the point of the initial static state (rz = 0,5, rm = 0.8) and the point of the new 
stable state (r= 0,877) we note that the time characteristic will pass through flow combination “p~— p" at 
first and then through “a— p",and the boundary value will be rz = 0.8 which corresponds to ry = 0.5. 


In the first section (*p— p") calculations are made according to formula (5), The integration constant 
is determined from the condition rm =0,8 and t=0, In the second section ("a— p") calculations should be 
made by formula (6) and the integration constant determined from condition rz=0,8 and t=, Here t; is the 
time of transition from the flow combination region “p— p” to that of “a— p", The calculation data is given 
in Table 1, 


TABLE 1 
*a-p" 
t, sec P, t, sec ra A, 
0 0.800 2.000 |t,= 6.72 0.800 1,250 
1.13 0.700 4.750 7.25 0.815 4.227 
2.76 0.600 4.500 7.84 0.830 1.205 
3.91 0.550 1.375 8.94 0.850 1.176 
t; = 6.72 0.500 1.250 12.25 0.877 4.140 


V = 462 cm? G = 980 cm/ sec” R = 2930 cm/degree 
T = 293° K, k=0.29, r=0.4, fe=0.004 cai? 


Figure 5, a gives the experimental curve of the transitional process, It also shows the calculated points 
(from Table 1), It will be seen from Fig, 5, a that the calculated points agree well with the experimental 
curve, 


Example 2, As distinct from the first example, the distance between the nozzle and the vane this time 
changes in the opposite direction, i,e,, jumping from hg = 0,083 to h,,,, =0.037 mm, It is required to deter- 
mine the time characteristic, 


In this example pressure Pyq in the interthrottle chamber at the initial static state will correspond to 
pressure P; .+,}, determined in the previous example, and the pressure at the new stable state will correspond 
to that of Py, also determined in the previous example, The rate of flow coefficient y of the second throttle 
should be taken at its maximum,* equal to 1 ,=0,48 (Fig, 3), Then, we shall obtaink=1, and f,=0,00113 


cm?, 

Py.kg/cm?® Py kg/cm? 
za} 
175+ 
ust 

a-p 
a | b 

u Pi, kg/cm? 

4} 145 

1 130 
a us 
c a-p P-P 
d 

Fig. 5 


Integration constants are determined from the condition that for “a— p” rg=0,877 at t=0 and for 
"p— p™ 1 =0,5 att 


All the calculation data is given in Table 2, 


TABLE 2 
*a-p* *p-p" 

t, sec atm t, sec P, = 2,5, atm 
0 0.877 1.140 = 1.36 0.500 1,250 
0.40 0.850 1.176 2.87 0.550 1.375 
0.78 0.830 41.205 4.85 0.600 1.500 
1.04 0.815 1.227 9.17 0.700 1.750 

t; = 1.36 0.800 1.250 18.00 0.800 2.000 


The calculated points are plctted in Fig, 5,b together with the corresponding experimental curve, The 
maximum deviation of a calculated point from the experimental curve is equal to 7%, 


* In view of the fact that yz is a function of the pressure in front of the nozzle, it would have been more 
correct in this case to have divided the whole range of P, variations into sections and taken the mean value 

of the rate of flow coefficient for each section, If, for the sake of simplicity, however, only one value of us 

is chosen, it should be the maximum (in this instance corresponding to the end of the transition process), since 
only in this case will the calculated value approach the experimental one, Physically this will mean that 
instead of taking the effective port area f, as increasing gradually during the transition period it is taken at 

its maximum value, It can therefore be expected that the calculated points at the beginning of the transitional 
process will lie below the experimental curve, 


ts 
e 
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Example 3, The distance between the nozzle and the vane is changed in one step from hg = to hgtab = 
=0, The feeding pressure is Py = 2 Ay, Pp = 1 Ag, T = 293° K, It is required to determine the time characteristic, 


It will be seen from above condition that in this instance it is the question of filling the dead-end 
chamber, Since r = 0,5 (Fig, 3) the flow into the chamber is pre-critical, The calculation should be made 
according to formula (11), 


Let us determine Py at the initial static state, The rate of flow coefficients of the first and second 
throttles can be assumed to be 0,8 since when fully opened the throttle represents a hole in a thin wall. 


Hence = dj/dj =0,045 and ry 1, =1 Ag, and =0,5, The area f,=0,00115 cm’, 
The integration constant is found from condition r, =0,5 at t=0. 


The calculated points have been plotted in Fig. 5,c and, as it will be seen from the graph, agree well 
with the experimental curve, 


Example 4, Let us determine the transient process in the interthrottle chamber of the relay (Fig. 3), 
when the feeding pressure changes in one step from Py =1 to 2,5 As, Pg=1 A, and h=0,05 mm, 


From the conditions of the problem it is clear that Py=1 A,, Let us determine the stable state pressure 
in the interthrottle chamber assuming that = 0.8. 


As a first approximation considering that Neb =1,2 A, we find from graph Fig, 4 that uh =0,.4 and 
= 0.737. 


According to the known ratio (f;/ f2)', considering that r = 0,4, we find from graph Fig, 2 that rp = 0.58 
and Ppl ap = 1.72 Ag. 


As a third approximation we obtain = 0.64, (f2/f2) =0.565, rz =0.68 and = 1,47 Ag. 


As the fourth approximation we obtain y iV =0.58, = 0.625, rp =0,.64 and Py staph = 1.56 Ag. 
The value of k is the one determined in the last approximation by ratio f;/f, = 0.625, 


Considering that we find f2= =0,00184 cm’, 


In order to establish the flow combination through the first and second throttles at the transient process 
it is necessary to take into account that at the point of the initial static state rp=1, r=1 and fy/f, =0, and 
in the new stable state r, = 0,64 and r=0,4, 


By following the curve r= 0,4, from the first point to the second it is seen that the initial stage of the 
transient process will occur in “a — p", the subsequent stage in *p— p™ and the boundary value of r, = 0.8, 
The calculation of the transitional process should therefore be carried out according to formulas (6) and (5), 


The integration constant for the case "a— p" is obtained from condition rz,=1 at t =0, and for the case 
"p— p” from condition ry =0,5 at t = ty. 


All the calculation results are given in Table 3, 


TABLE 3 


0 1.000 1.000 1. 

0.52 0.950 1.052 7.22 0.550 1.375 

1.24 0.900 1.110 14.81 0.600 1.500 

1.98 0,850 1.180 . 24,23 0.6.0 1.550 

3.02 0.820 1.220 33.03 0.625 1.562 
0,800 1.250 


Points found by calculation are plotted in Fig. 5,d, The same graph for comparison carries the experi- 
mental curve, 


a t,sec | Tr | P, = 1 | Tu As 
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EMERGENCY CONNECTION OF LAMPS 


B. S. Sotskov 


(Moscow) 


A method of emergency connection of lamps (electric supply, lighting, etc.) by 
placing them in parallel is examined and the required values of the voltage and ballast 
resistance are determined for the condition that when one of the lamps fails the re- 

sulting optical and emission characteristics remain unchanged, 


Many automatic devices use incandescent lamps and electron tubes, The electron tubes are used as 
amplifiers, oscillators, detectors, relays, etc,; incandescent lamps are used as sources of light for photoelectric 
amplifiers, photorelays, photoindicators, etc, 


Since the lamp filament under the action of the heating current wears and burns out in time, emergency 
connection of lamps become necessary, The simplest method of emergency connection is the parallelling 
of two lamps, A simple parallel connection however, is not suitable, since the failure of one filament changes 
the total anode current in the case of tubes, Similarly in case of incandescent lamps 
the resulting luminous flux is changed, Above circumstances lead to a search for other 


Go methods of emergency connection, 
by A constant anode current or luminous flux when one of the parallel filaments 
fails in the case of electron tubes and incandescent lamps, respectively, can be 
assured by means of the circuit connection shown in Fig, 1. 
Let us examine the connection of two lighting or signalling lamps with the view 


Fig. 1 of increasing the reliability of their work, Let there be a fixed condition that the 
lighting of the surface must remain constant whether one of two lamps are alight, 


Let us attempt to find the required condition of operation for the lamps, 


Figure 2 shows the relation between the luminous flux and voltage across a lamp @ = f (U) [more 
accurately / 19) = f (U/Uyq)] in the form of curve 6, The same graph shows the relation between the 
current I through the filament and the voltage across the lamp [more accurately I/Iy, = f (U/Ujo9)] in the 
form of curve I, The graph also includes curves for the luminous flux and current for two lamps (curves 26 and 
2), Here Ugg and 499 denote the nominal values of voltage and luminous flux for a lamp, 


Let us assume that a luminous flux @, is required from the lamps, Let us find on the curve 24 point a, 
corresponding to the required flux , (or x/499). It corresponds to voltage Ug (or U,/Ujg9). The same 
voltage U, corresponds to a current through the lamps determined by point a‘ on curve 21, 


If one of the lamps burns out, the luminous flux according to the operation condition, must remain con- 
stant, For the value $, we obtain on curve @ a point b, It corresponds to voltage Uj), (or U}/Ujg9) at which 
the current through the lamp is determined by point b’ on curve I, 


In order to obtain these relationships the lamps should be connected as shown in Fig, 1, In order to find 
the value of resistor ry and that of the supply voltage E, let us draw a straight line through point a and b, This 
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line will cross the X-axis at point c which determines the required voltage E, (or Ex/Ujgg), Resistance to is 
U 
determined from condition ae am =t (where Uj} and I, correspond to point b* on curve I), 

In the case of the operation of lamps with photocells 
it is necessary to know the relationshipI,= ¢,(U), where 
I, is the current through the photocell at a given position 
of the lamp and the photocell, and U is the voltage at the 
lamp terminals, By using relationship I, = ¢ ;(U) instead 
of the curve 6 = f (U) we obtain in a similar manner the 
values of rm and Ey (Fig, 3), Curve 21, for the photocell 
current when it is lighted by two lamps can be found in 
various ways, Figure 4 shows the luminous characteristics 
of photocells with a constant color content of the luminous 
flux and with the voltage at the photocell = const, In 
the case of a linear light characteristic of a photocell 
1, = ¢ (U) (curve a, Fig, 4) it is possible to double the 
ordinates of the curve I, = 1 (U) which will produce the 
required curve 2 = g,(U), In the case of a nonlinear light 
characteristic of the photocell (curve b, Fig. 4), it will be 
necessary to determine curve 26 = ,(U) experimentally 
with a luminous flux supplied to the cell from two lamps, 


2/_=9,(U) when lit by 
/ two lamps 


view 


he 
> and 


» find 


>. This Fig. 4 Fig. 5 
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Instead of obtaining the relationship I, = y ,(U) and 2lg = gy 2(U) experimentally, it is possible to 
find them by calculations, 


The above technique of calculating the supply voltage and ballast resistances is also applicable to 
parallel connection of electron tubes, 


In their case it is necessary to plot the relation 

Zt between the anode current and the heating voltage I, = 
= F,(U,,;) at U, =const and U, =const, where U, is the 
anode voltage and U,, the grid voltage (Fig. 5), and the 
relation between the heating current and the heating voltage 
= f, (Up). By plotting curves 2g =F,(Uy) and = 
= f,(U}) for two tubes connected in parallel and curves 
I, and f,(Up) for one tube (Fig, 6), we 
obtain for a given value of anode current points a and b 
on curves I, =F,(U,,;) and 2, =F and the values 
of voltages Ujjg and Uy. Having next found points a’ 
and b’ for voltages Uzj, and on curves and 
a sm line is drawn through these points,and E, and 

x_— 
sponding to point b*, 


ro= are found; here Ip; is the current corre- 


Depending on the matching of the characteristics 1, = F,(Uyy) and I}; = f (Uj) of each of the actual 
tubes with the corresponding calculated characteristics,the accuracy in maintaining a constant Ig will be 
within the limits of 1,5 to 8%, 


The lifetime of parallel connected tubes increases since, in order to obtain the same anode current (or 
luminous flux in the case of incandescent lamps) a smaller voltage is required. The relation of lifetime T, 
to voltage U, can be determined from formula T, = Ty (Uyo9/Ux)% where a ss 8 to 12,depending on the 
construction of the filament, and Ty is the lifetime atthe nominal voltage Uj. 


The combined lifetime of two lamps connected in parallel is equal to Tys =Tyg+T,p, where Tyg 
and T,, are the lifetimes at voltages U, and Uj, or 


U a U 2 
+ 
The value of a for each tube is not the same and is usually subject to the normal distribution law, 


Usually the mean value is taken for a, In some cases it is expedient to determine two lifetimes using the 
values Crpin 29d Gmaxe This will provide an evaluation of the possible spread in the lifetime, 
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SIMPLIFIED ALGEBRAIC SYNTHESIS OF RELAY CIRCUITS 


la, I. Mekler 


(Moscow) 


The paper deals with a simplified algebraic synthesis of relay circuits according to 
a method suitable for obtaining a structural formula directly from switching table, thus 
avoiding all the operations and subsequent simplifying transformations employed in normal 
circuit design methods, A comparison is made between this and other known methods, 


The formulation of operating conditions for various elements of multiphase circuits and the compiling of 
algebraic expressions which determine these conditions is accompanied by certain difficulties due to the 
necessity of analyzing at the same time the internal interactions in the circuits. 


It is known from the theory of relay- contact circuits [1] that compiling algebraic expressions which 
characterize the operating conditions of a circuit and its elements can be carried out by the following methods, 


1, By copying out from the switching table all the combinations of element conditions at which the 
given element must be in a connected (disconnected) condition, This is explained by the fact that each com- 
bination of element conditions for which the circuit of any element must be closed (open) is the condition of 
its operating (not operating).* 


2. By means of copying from the switching table the condition of elements in the phase preceding the 
one in which the given element is connected and in the phase preceding the disconnecting of this element, 
Since the condition of elements in the phase preceding the disconnecting of this element is the condition for 
opening of the element's circuit the algebraic expression of this phase must be taken with an inversion. 


The algebraic expressions compiled in the usual manner become complex and contain as a rule a large 
number of duplicated and mutually excluding circuits. It is true that the technique of the relay- contact cir- 
cuit theory permits one to reduce the initial complex algebraic expressions to simpler ones, but this is 
panied by a complicated process of transforming the formulas originally obtained, 


The method [2], worked out by M, A, Gavrilov, of dividing one unwieldy switching table into several 
simpler ones only consisting of elements which react on the element in question, rather simplifies the algebraic 
representation of the functioning of circuit elements but at the same time lengthens to a certain degree the 
synthesizing process. 


Even when the above operations are applied the problem of compiling 4 circuit remains fairly complicated 
and requires further transformations, 


The method [3], proposed by A, N, lurasov, of obtaining a structural formula from the switching table 
requires subsequent complicated checking and has not, therefore, been gener2ily adopted, 


* In this instance compiling algebraic expressions for circuits based on conditions of nonfunctioning is in view. 


4 
+ 
. 
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The method described in this article is based on the fact that the elements comprising the circuits have 
a different effect on any one element Xj, depending on the position these elements occupy in the switching 
table, By examining the conditions under which these elements or groups of elements affect element Xj it is 
possible to establish the rule according to which these effects are produced, The rule of interaction between 
elements differently placed in the switching table thus obtained permits one to avoid complicated transforma- 
tions of structural formulas and to obtain algebraic expressions in the final or nearly final form directly from 
this table. 


1. Schemes Requiring Discrete Circuits to Retain Blocking Arrangements 


Through Their Own Contacts 


The functioning condition of element X; operating once per cycle is written as follows [1]: 
Ix; = F,+ (1) 


where F, is the Boolian function representing the product of variables which correspond to contacts of circuit 
elements (with the exception of element X;) in the phase preceding the functioning of element Xj, F2 is the 

Boolian function representing the product of variables which correspond to contacts of circuit elements (with 

the exception of element X;) in the phase preceding the releasing of element Xj. 


Since the discrete elements of the circuit controll- 
ing element Xj, under consideration, occupy different 
annem --------- positions in the switching table with respect to above ele- 


Connected condition of element Xj 


JJaversion._ - ment and hence affect its operation in different ways, it 
TA) lesen oe is expedient to divide these elements into categories with 
respect to the functions they fulfill (Fig, 1).* 
[A.J Ag is the aggregate of circuit elements which change 
re A; 7 their condition once during the interval between the 1-st 
and 2-nd changes in the condition of element Xj;, under 
inversion [ay Jlnverdion 


Ag is the aggregate of circuit elements which are 
switched in and out before the 1-st change in the condition 
of element Xj, or after the 2nd change in the condition of this element; | 


Fig. 1 


Ag is the aggregate of circuit elements which are switched in before the 1st change in the condition of 
element Xj and switched out after the 2nd change in the condition of the same element, 


Ag is the aggregate of circuit elements which change their condition twice during the interval between 
the 1st and 2nd changes in the condition of element Xj, under investigation, 


The number of elements in categories As3,A 4, As, and Ag is arbitrary. 


Category As elements. If side by side with elements of other categories there is but one element of 
category Ag,the finding of the position of its contacts does not present any difficulties. In the cases when the 


number of elements of this category is greater than one, and in practice can be any number depending on the dl 
problem in question, it is necessary to take into consideration the order of their switching in (Fig. 2). 
Let us denote by fs a Boolian function representing the product of variables which correspond to the 
contacts of elements of the Ag category in the Ist condition 
X 
where Xg.4Xg.2. » . Xge are variables corresponding to contacts of category As elements in the 1st condition 
* In order to avoid indexing errors the numbering of A subscripts follows those of F, fo 
th 


1100 


4 


and placed in the switching table in the order of increasing subscripts*; fs is a Boolian function representing 
a product of variables which correspond to contacts of category Ags elements in the 2nd condition 


— 


Ie HE (3) 
It should be noted that f*s is a Boolian function re- 
| Connected condition | presenting the product of variables which correspond to 
i___of element Xj ls contacts of category Ag elements in the phase preceding 
the releasing of Xj. 


' a) The condition of contacts in the retaining circuits 
will be expressed by an inversion of fs, i.e., 


Connected condition of 
€lement X; circuit | = + +... + (3°) 
Fig. 2 In the most general case elements of the Ag category 
a are connected (disconnected) in the order of increasing 
subscripts and disconnected (connected) in the same order, 


Hence, it would appear, that the algebraic expressions of circuits consisting of these elements are not 
simplified by the formulas for elements functioning in a definite sequence, 


If, however, the functioning of these elements is considered in conjunction with X;, expression xjXg.4 + 
+ XjXg.2+.. .+X;Xgen is seen to represent a certain sequence in the functioning of these circuits, namely: 
simultaneous connecting, but disconnecting in the order of increasing subscripts, 


Hence, 
+ +... + = (4) 


Thus, in examining circuits with more than one element of the A; category it is found that in retaining 


circuits there remains only one contact x3. which is the last one in the order of element distribution in the 
switching table, 


Category Ay elements, The Boolian function representing the product of variables which correspond to 
contacts of category Ay elements in the phase preceding the releasing phase will take the form 


hs = (5) 


In this instance these elements can be omitted from the algebraic expression of the retaining circuit, 
since X; and Ag are not connected simultaneously and therefore 


+ +--+ + = 0. (5") 


Category Ag elements, The Boolian function representing the product of variables which correspond to 
the contacts of category Ag elements in the phase preceding the releasing will take the form 


In this case these elements can be omitted from the algebraic expression of the retaining circuit, since 
X; and Ag are switched in the reversed and switched out in the direct order, and hence 


+ + = O. (6") 


* Since in the most general case the circuit contacts can be either closing or opening, their symbols in the 
formulas must be accompanied by the sign ~, denoting that the contact can be either one or the other, For 
the sake of simplifying notation this sign has been omitted, 
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Category Ag elements. If the number of category Ag elements is larger than one they can occupy in 
the switching table different positions with respect to element X; and to each other, for instance, in the manner 
shown in Fig, 3, 


Fig. 3 


Let us denote: X¢.3, Xge2, - . - » Xgen to be variables corresponding to contacts of category Ag elements, 
at the 1st change in condition of Xj, placed in the switching table in the order of increasing subscripts; Ts 
to be a Boolian function representing variables corresponding to contacts of category Ag elements at the Ist 
change in condition of Xj: 


he eee (7) 
In the algebraic expression of the phase preceding the releasing, the symbols of contacts of the category 
under consideration are represented as multiples, i.e., the contacts in the circuit are connected in parallel. 
In examining Fig. 3, a it will be seen that 


Te = + + Zeon = (7") 


according to the principle of the equipollence of elements functioning with a given sequence, viz, connecting 
in the order of subscripts and disconnecting in the reversed order, 


If the elements take up the position shown in Fig, 3, b where they are connected and disconnected in 
the order of subscripts, the algebraic expression cannot be simplified in the normal way. Below we give a 
method of simplifying a circuit with such a sequence of operation of its elements. 


T 


Fig. 4 


It should be noted that if in a circuit controlling some element Xj there are only category Ag elements 
present, this circuit cannot be set up, in view of the presence of contradictory conditions: the circuit elements 
in this case must have the same conditions in the phases where X; is connected (at the beginning of connecting) 
and in the phase preceding the releasing, as well as in the phases of a released condition, This can be easily 
seen in Fig, 3. 


In order to make the circuit applicable it is necessary to introduce an element of category As. The 
presence in the circuit of an element of this category is a necessary condition of its applicability. * 


* An element of category Ag is a necessary condition for setting up the retaining circuit, Since the realiza- 
tion of the scheme also depends on the ability to set up the connecting circuit, the above condition is necessary 
but not sufficient, 
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After inserting in the scheme an element of the category Ag (x3.n) the graph under consideration will 
assume the form shown in Fig, 4, a, Here the elements X¢.1, Xg.2, . . . Xgen can be considered as being connected 
simultaneously and disconnected in the order of their decreasing subscripts, 


The algebraic expression of the retaining circuits of such a scheme can be simplified: 


The presence of a category Ag element in the switching table, as well as making it applicable, solves 
the problem of simplifying the scheme of Fig. 3, b, In this case, as it will be seen in Fig, 4, b, elements 
» Xgen CAN be considered with the presence of as being connected simultaneously and dis- 
connected in the order of their increasing subscripts, The algebraic expression for a retaining circuit of such 
a scheme will be 


+ + + Lon) = Li + Leen): (8") 


So far only mutual positions of category As and Ag elements were considered in which the significant 
element of the Ag category (its contact Xg.n) fell within the zone occupied by all the category Ag elements, 
It is of interest to examine the case when the 
- 1X; element of the Ag category falls in the zone occupied 
+ 'Z5y by only a part of the category Ag elements, For in- 
———“ 25.2 stance, let us examine the graph in which contact 
we nee Xgen remains connected after the switching in of 
contact and is disconnected before the switching 
in of contact X¢.2 (Fig. 5). 


‘730 In this case in a manner similar to the pre- 
Fig. 5 ceding one, the algebraic expression for the retain- 
ing circuit will be 


Li + + Loon + + Leen) = Li + + Zen). (8") 


It should be noted that the circuit elements can pass in the course of operation from one category to 
another, For instance by examining Fig. 6, a and b, which give the switching table for a consecutive opera- 
tion of the circuit under the effect of multiple incoming pulses, it will be seen that the element x of category 
Ags connects Xj and then passes to category Ag. Figure 6, a shows a sequence of circuit operations in which 
the condition of contact x is the same in the phases preceding the switching in and out of Xj, whereas in 
Fig. 6, b the condition of contact x is different in the two phases, Lines on the diagram of Fig. 6 indicate a 
functioning condition of the appropriate circuit elements. 


An algebraic expression of the circuit operating conditions with respect to both cases is given below, 


Let us denote the contacts of element x when they comply with conditions of category Ag by x3., and 
when they comply with those of category Ag by Xg,,,. 


| Connected condition | | Connected condition | 

of element X i | of element X 

en 

| 

| Connected condition of | Connected condition saad 
| the circuit | the circuit 


a b 


3 

| 

Fig. 6 
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According to above reasoning Expression (1) can be replaced by expression 


(9) 
+ (Xn + + + Len + 
It is known from [1] that Expression (1) can be given another form 
Ix, = (Fi Fo. (10) 
This form can be represented in another way 
Ix, = + + F + , + 
+ + Lor + +--+ + + = 
= (Fy + 2%) + + +... + + (11) 


In Expression (11) F F =0. If = then =0. Also F = F;. Let 


Finally, we obtain 
F = (Fy + 2) + 20-1 + 2-2 +--+ + %,.,)- (12) 


In comparing Expressions (11) and (12) it will be seen that the number of symbols of elements in the 
latter is smaller than in the former. 

Hence of the two versions X3.; = Xba. and X3.1 = Xg,,_ the former is preferable for compiling switching 
tables, 


2. Schemes in Which Contact Combinations in Phases With A Connected 


Circuit State Are Not Repeated in Phases With A Disconnected State 


In deducing Formulas (9), (11), and (12) we considered the most general case according to which the 
setting up of the scheme in the retaining circuit with respect to any element required blocking by means of 
its own contact x;. 


Let us denote the sum of element conditions, in phases with a connected state of circuits, affecting 
element X; by ZF let us separate from these circuits the contacts of this element, as it was done in [2], and 
denote the remaining expression by LF3, 


Let us examine a case when not a single combination of element contacts contained in LF§ i.e., 
corresponding to the connected state of circuits affecting element Xj, is repeated in the sum consisting of cir- 
cuit-element contact-combinations corresponding to the phases of a disconnected state of circuits of the same 
element Xj, a sum we shall denote by =F p, and after separation from these circuits of the contacts of element 
X; denote it by <XF',. With such a relation of element conditions in phases of a connected and disconnected 
circuit the conditions of functioning of an element controlled by them can be written in the form 


tx, = (13) 


It will be seen that the requirement of blocking through its own contact x; does not exist here. 


According to Expression (13) the condition of operation for element X; is represented by the sum of the 
circuit element conditions in all the phases when X; is functioning or, which is the same thing, by the sum of 
unit constituent elements in the phases when the above-mentioned element X; is affected by the functioning 
of other elements. 
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Let us denote this sum by rk. 
In this case Expression (13) can be written in the form [5]: 


fx, = UF; = Kf. (13") 


In a manner similar to the one described above let us examine circuits DF‘, which consist of various 
categories of circuit elements depending on their position in the switching table with respect to the element 
X, under consideration, Elements of categories Ag, Ag, As, and Ag should be examined, 


Category Ag, Ag, and Ag elements, If there is only one category As element in the circuit there is no 
difficulty in finding the position of its contacts in the circuit, If, however, in addition to elements of other 
categories there are in the circuit several elements of this category and their number can reach any value, it 
becomes necessary in solving this problem to take into consideration the order of their disconnecting (operation), 


The symbols of category As element contacts are contained in EF3 in the form of a sum, and the 
corresponding Boolian function according to the notations adopted above can be written as: Xg.4+Xg.2+...+ 
+Xgen= F's. 


In the most general case it is necessary to examine elements of this category which are connected (dis- 
connected) in the order of their increasing subscripts and disconnected (connected) in the same order, and 


which, it would appear, are not subject to simplification by means of formulas for elements functioning with 
a definite sequence, 


If, however, one separates from the circuit under consideration element Xq which changes its condition 
in the phase preceding the one in which element X; changes to its functioning condition, and then examines 


the operation of As elements in conjunction with element Xq, one will observe a definite sequence in the 
operation of the circuits, 


It should be noted that element Xq can belong to the Ag, Ag, OF Ag categories, 


Connected condition | |. 


of element X; of element X 
1, _ : 
+ 

‘ 

| Connected condition of i, — 

the circuit Connected condition of 
| | | the circuit 
a b | 


Fig. 7 


If elements of categories Ag and Ag are taken as Xj it is sufficient for the solution of the required pro- 
blem to represent Expression (13) in the form 


where rk) is the sum of unit constituents of the phase by phase condition of circuit elements, in which xq 
is missing, Ki) the zero constituent equal to EKp}?, 


* The author has proved that such a sum of unit constituents is equal to the zero constituent of the elements 


which change their condition during the functioning condition phases of element X; in question, Let us de- 
note this constituent by ), 


|| 
fx, = = EK} = a4) 
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Inspection of the graph in Fig. 7, a shows that the category Ax elements considered in conjunction with 
Xg Obey a law of a definite sequence, namely simultaneous connections, but disconnections in the order of in- 
creasing subscripts 


The full condition of functioning for element X; will be according to (14): 


Ix; = (Z3.n 
Figure 7,b shows simultaneously and .4, Xs.2, . » Xgen- 
Similarly to the above reasoning +. . Xg.nXq=XgenXq and f xX; = Xq (Xg.n + 


Thus, it will be seen by inspection of circuits with category As elements which do not change their 
condition during the phase preceding the change of condition in X; that the algebraic expression of the circuit 
operating conditions contains only one contact Xs.n, which is the last in the order of contacts in the switching 
table, 

If elements changing their condition in the phase preceding the change of condition in element X; 
should include category As and not category A, or As elements, then only one of the category As elements 
will be included in the conditions of operation for element X;, There will be no other elements of this cate- 
gory in the algebraic expression for the circuit, 


This can be easily seen in Fig. 8, a which shows that the algebraic expression for the circuit is fx, = 
= + Xe. 


| Connected condition of | | Connected condition of 
| element X; | element xj 
+ | 
— 25 n | 
} Connected condition of “ Connected nditi f 
the circuit the cireutt 
a b 


Fig. 8 


According to Fig, 8, b the algebraic expression for the circuit can include only X3.; and xg, Other 
elements of the category Ag such as X3.2,.... , Xgen Will not be included, Moreover, the circuit cannot 
be set up at all without a blocking action by means of one of its own contacts. 


The relation of category Ag elements with elements of other categories has been discussed before, The 
previous statements with respect to them hold good in this case as well and do not require any additional 
analysis, 


It was previously pointed out that circuit elements can change categories in the course of their operation, 
Figure 6 shows an element of category As which switches in X; (it is denoted as xs.) and then changes to 
category Ag and is denoted as X65.4° Above phenomenon can occur in this case as well. 


On the basis of above statements the algebraic expression of the conditions of operation for element Xj, 
written above in the form of a sum of contact conditions of elements affecting X; in all the phases of its 
functioning state, will be as follows. 


1, If an element of the jy or As category changes its condition in the phase preceding the functioning 
phase of X; under consideration: 
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2, If an element of the Ag category changes its condition in the phase preceding the functioning phase 
of X; under consideration: 


Ix, = %s + + + (15') 


Here xs has no subscript since it is a contact of the As category element which changes its condition in 
the phase preceding the functioning phase of X; under consideration, 


In deducing Formula (14), and in determining the number ot element contacts of the As category which 
are included in the algebraic expression of the operating conditions of X;, only one element of the A, or Ag 
category was examined, This corresponds to the conditions in a large number of applied circuits, In the most 
general case, however, the number of elements of the said categories can be increased indefinitely, Since 
these elements do not change their condition during entire period of the functioning state of X; they can be 
taken outside the summation sign of unity constituents and become factors of the zero constituent: 


where K\) is the sum of unity constituents of the phase-by-phase condition of circuit elements which lack 
all the elements of the A, and Ag categories; and K@) is the zero constituent equal to 5K) . 


The algebraic expression of the functioning conditions of element Xj can be written for the most general 
case in the following form: 


fag = (9.0 + + + + 
Moreover, if,in the phase preceding the functioning phase of X; under consideration an element of the 


Ag OF Ag category changes its condition, then xs.n is the element contact of the Ags category, which is the last 
in the order of their position in the switching table, 


If, however, one of the elements of category As changes its condition then x3.p will be the contact of 
that very element, 


3. Examples 


Examples of setting up operating conditions for all circuit elements by means of an algebraic method 
according to Formulas (9), (11), and (12) and according to other methods are given below for comparison pur- 


poses, 

Phas 0 i 2 3 4 5 6 7 8 9 |} 10 | if 12 | 13 | 144 | 15 | 16 
—ai+a —a +a —a +a —a +a —a 

—Xi} 

m 
+X5 
Fig. 9 


Let us inspect the switching table given in Fig. 9 and determine conditions of operation of circuit elements 
by former methods described in [1] and [4] and by the proposed method according to Formulas (9) and (12). 


4 
. 
y 
= 


Above examples of comparison between circuit synthesis by the ordinary methods and the method pro- 
posed in this paper show the advantage of the latter, This method does not require complicated transformations 
and provides results directly from the switching table, 


Methods according to Formulas (9) and (12) 


Methods according to[1 and 4] 


+ + + +24) 
Since X, and Xz as well as X; and X, are not 
connected simultaneously 
=0 = 0 
+ + 29) 
Or 
{X,== + + x(a+ = 
= + - + 22) 
2) ax, + = 
== + + 2, + + 
Since ¢ 22% + + = 
(circuits connected with a definite sequence) 
= + (a + 24) 
or 
+ + (a + 24) = 
== + (a + 2%) 
3) fx,= + = 
== 12,24 (a + + 24) 
Since and are not connected 
simultaneously = 0, 
+ (a + + 24) 


fx, = + ax + ax 
+ + + 2%) = + + 23) (a + x2 + %) 
== AX, + + 2, + + Zs) 


Since X, and Xy are not connected simul- 
tancously xx, = 0. 


Since Xs and X, are connected in the order of 


increasing subscripts but are disconnected in 
the reverse order 


1) a, fa = 2s, f6,., 
According to (9) 


1x, = (axgxq + 2) (a+ 


2) Is = 4, v1, Zs, fos 
According to (9) 
Gx, + 22 (a + 2%) 
According to (1.2) 
+ 29) (a + 2%) 


1, 
According to ( On 


{X,= + (a + + 
According to (1 “4, 


fx,= + (@ + + 2) 


4) a, 22; fa= 
fan. = a 
According to (9) 
OT +24 (a + 2) 
A 


_ = 0 
és fy, == + (a + 
+ ax,x3a + 22) + 
+ (axyxs + 24) (a + 23) 


It has been stated previously that in addition to the above method there are other methods which do not 
require copying the structural formulas for all the elements contained in the switching table, but only require 
determining the elements which are in the circuits of the element whose conditions of operation are being 
determined, 


Among such methods is the one proposed by M, A, Gavrilov [2] which consists in dividing the basic 
switching tables into auxiliary ones, 


The previously used examples are now solved by means of compiling new switching tables, only con- 
taining elements which change their positions in the phases preceding the functioning or releasing of the 
element under consideration and then checking the applicability of the tables thus obtained, 


The new switching tables are given in Figs, 10, a ; 11, a; 12, a and 13, a, It will be seen that not one 
of the newly compiled tables is applicable, and according to the rules of the method it is now necessary to 
find elements in the general table which would make the tables applicable, and add these elements to each 
new table, 
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B 2° |-al+a —a +a —a +a —a +a —a 
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I 2° |—-al+a —a +a by +a —a +a —a 
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22 |—X, 
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c 
Fig. 10 
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» sofa |e |e | 42 15 | 46 
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a 
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Such elements were found to be the following: X, and Xg in obtaining operating conditions for element 


Xy; and X, for element Xz; X,, Xz, and X, for element and X_ and Xz for element X,. 


in Fig. 12,d and for element X, in Fig. 13, c. 
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Thus, the operating conditions of the circuit elements are determined as follows: 


The determining of the additional contacts does not present any difficulties and is shown in Figs, 10-13, 
The final switching table for element X, is given in Fig. 10, c, for element X, in Fig. 11, c; for element Xs 


fy, = + = + 21 (a+ 29), 


fy, = + + (a + 2%), 


= + == + (a + + 2), 
fy, = = = + 24 (a + 
It is useful to inspect the switching table given in [2] and to compile operating conditions for circuit 


elements by the method recommended here according to Formulas (9) and (15) and to compare it with M, A, 
Gavrilov's method. 


The full switching table is given in Fig. 14, 


5 iat? 9 | 10 

2° |—a\+a 
Binary, XxX, +X, +X, 

2? +Xe | 

23 

Sum of ay y ew 
combina- | 0/1/ 3 718th Pai 3 

tions 
Fig. 14 


It is shown in [2] that according to M. A. Gavrilov's method in order to determine the elements affecting 
respectively X,, Xp, Xs, and X,, it is necessary to compile tables which form parts of the table in Fig. 14. It 
is also stated there that in determining operating conditions for each element these tables can turn out to be 
invalid, as is the case in determining the conditions of operation for X;. This fact leads to the re uirement of 
inserting an additional element, that is, to compiling a new table, The determination of the operating con- 
ditions for elements X,, Xq, Xs, and X, by this method is given in [2]. 


Let us determine the conditions of operation of the same elements X4, Xg, Xs, and X by the proposed 
method according to Formulas (9) and (15). 


The conditions of operation for element are determined by Formula (15): fs = 
FX, =a xg). 


The conditions of operation for element X, are: fy=%q, %, fg =Xs. According to Formula (15) fX,= 
= 


The conditions of operation for element Xs are: f3=%, fg =X. According to Formula (15) = 


The conditions of operation for element X, are: fg=Xg, X,, Xg. According to Formula (9) fX,= 
= AXy t+ Xp 


Taking into consideration formulas for a given sequence applied to connecting circuits (the retaining 
circuits obtained by Formula (9) are not subject to simplification) we obtain fx, = Xs + X Xp. 
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SUMMARY 


1, The method described in the paper lends itself to a rapid calculation of algebraic expressions for 
circuits direct from a switching table without all the operations and the accompanying simplifying transforma- 
tions which are employed in the existing methods of setting up circuits. 


2, M, A. Gavrilov's method, based on the division of the switching table into several simpler ones with- 
out changing the principle of setting up the circuit, simplifies the synthesizing process but at the same time 
rather lengthens it, 


3. In comparing M, A, Gavrilov's method with the one described in this paper, it should be noted that 
the determination of the connecting circuits is accomplished by the former quicker than by the latter, It is 
recommended, however, to carry out the determination of the retaining circuits by the method described in 


this paper, 


LITERATURE CITED 


[1] M, A. Gavrilov, The Theory of Relay Contact Circuits (AN SSSR Press, 1950). 


(2) M. A, Gavrilov, “Basic formulas for synthesizing relay circuits," Avtomatika i Telemekhanika 
15,6 (1954). 


(3) A. N, Iurasov, Setting up Structural Formulas for Multiphase Circuits [In Russian] Moscow 
Mechanical Institute, (Mashgiz, M,, 1952). 


[4) Ia, 1, Makler, Constituents Sufficient for Ensuring the Operation of Relay-Contact Circuits (Trans, 
AN SSSR 117, 4, 1957). 


Received June 8, 1958 


1114 


. 


AN OPERATIONAL AMPLIFIER WITH A DIFFERENTIAL INPUT 


Vv. B. Smoloy 


(Leningrad) 


The paper describes block schematics of amplifiers capable of simulating, 
due to specially designed input and feedback circuits, the greater part of integro- 
differential relationships important from the practical point of view, 


Operational dc amplifiers with a compensation circuit (see sketch) can be used for solving electrical 
simulation problems. 


Operational amplifiers operating on this principle were 
proposed by Philbrick in 1949 and form the basic component 
6s, of an electronic analog made by a U, S, firm of the same 
a : as name, These amplifiers are capable of simulating, due to 
% the differential input circuit, a wide range of functions, 
bt, Some important practical applications of these amplifiers 


4,” cc ra produce simpler circuit designs than conventional Soviet 
made amplifiers used in electronic simulating devices, 


In order to demonstrate the possibilities of an opera- 
by, 8 by tional amplifier with a differential input let us determine 
7 G the relationship between its output voltage U,,, the input 

voltage and the conductances and 


vm From inspection of the block schematic it immediately 
Y follows that the output voltage Uw is the root of equation 


U 


A substitution of junction voltages Ua and Up by input voltages Ux, and Uy, and the output voltage 
determines the (mplicit function 


i 1 


G G 
(p) = hy, (P) = 
D (P) Gy (P) + P) 
1 


m 
and transfer functions of the corresponding input circuits, and h,, (p) = 1—'h,, (p) 
i 


a 
Uy, 
where 
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TABLE 
Amplifier block Simufat latignshi 
Amplifier operstes 
m 4 : summator 
1+R 
R 
n 1 
k 
py Rex 

R 


on” xk — 


Amplifier operates 
as an algebraic 
summator and 
differentiator 


Amplifier operates 
as an algebraic 
summator and 
integrator 


Re, 


Amplifier operates 
as an algebraic 
summator, 


differentiator and 
integrator 


= xx + 
n 
Uy, 1 
xk 
Ly 1 
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By = Cyx, Ax = Aix 
1 
n 
— 1 
a +p" b Ag U |, 
Y, Ay mn 4 
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n m 
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Cokin 


TABLE (Continued) 


Amplifier block 
schematic 


Simulated relationship 
AU xk, Uyk, (P), hyx(P)] 


Remarks 


n m 
1 1 
ak 


= k 
1 


—U pt, hep = 
v>) okP + "xk 
1 


y= 
AaxP 
1 
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m 
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1 


h 


Algebraic summator 


and differentiator 
of the m-th order 
bye and by 
are transfer functions 
of the input and 
feedback circuits) 


n 
+ 
1 


—k m = 
1 1 


1 


n 
m—1 


i+ B,P* 
1 


(k=1, 2,--+5 n), 


Algebraic times m 
summator and 
integrator (hy 

and h,, are 
transfer functions of 
the input and feedback 
circuits) 


Algebraic summator and 


integrators with U, = 
U.. =U the simulated 


relationship is U,, = 
=U 
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is the transfer function of the feedback circuit, 


When condition | why(p) | >» 1 is fulfilled the final expression for the amplifier output voltage will 
take the form 


n 
U. = 
(P) 


It foliows from the last expression that for input voltages Uy, the amplifier acts as a repeater and does 
not improve the quality of the passive summation network U,, but for voltages Uy, the amplifier represents 
a normal anode to grid feedback circuit and considerably improves the quality of the summation network Up, 


This circumstance must be taken into account when the amplifier is used in electronic simulating de- 
vices, 


The shape of the simulating characteristic of Uy is determined by the type of conductances G(p) and 
the way they are connected in the input and output circuits of the amplifier, 


The table shows the most typical practical schematic circuits of this amplifier for simulating automatic 
control systems we well as relationships corresponding to these circuits, 


The accuracy of simulating relationships given in the table is determined not only by the theoretical 
error 


G 
100 2 


but also by the instrument error, due to the difference in the amplification factors 1 and pi of the first 
differential stage of the amplifier. This error is given by the formula 


i Avy 


where Ay, and Ay, are quantities representing the instability of the amplification factors of the first stage 
tubes, 


At 


= 


3...) 


1 1 


SUMMARY 


1. An operational amplifier with a differential input is a universal device capable of simulating re- 
lationships (algebraic summation, integration, differentiation and other more complex integro- differential 
functions) without changing the sign of the input voltage, As a result of this property it becomes possible to 
omit in analog circuits special sign- changing (inversion) operations which require additional amplifiers with a 
1, 


n 
(P) 
1119 


2, The accuracy of the described amplifier is a little lower than that of conventional operational am- 
plifiers due to the error in obtaining the voltage difference in the differential stage of the amplifier. However, 
the schematic simulation of the majority of function given in the table by means of conventional amplifiers 
is more complicated, Hence, the possible decrease in accuracy when a differential input amplifier is used 
should be considered on the merits of each particular case, 


8, With zero input voltages U,;, = 0 the amplifier operates as a conventional anode to grid feedback 
amplifier, 
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FIRST IFAC CONGRESS 


The International Federation on Automatic Control (IFAC) will hold its 


first Congress in Moscow in 1960, 


y The USSR National Committee on Automatic Control (Moscow, 15 a Kalanchevskaia Street) invites all 
| specialists engaged in scientific and practical work in the sphere of automation to prepare for this Congress, 


Persons wishing to read papers at the Congress should apply to the National Committee, enclosing the 
title of the paper with a brief annotation, before February 15, 1959. 


The scientific program of the Congress and its organizational details will be published in the journal, 
Avtomatika i Telemekhanika No, 1, 1959, 


USSR National Committee on 
Automatic Control 
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